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CHAPTER  I .  INTRODUCTION 


► 


1.1  Magnetostatic-Wave  (MS W)  Devices 

The  purpose  of  this  work  is  to  study  magnetostatic-wave 
propagation  in  a  rectangular  waveguide  partially  filled  with  a 
lowloss  ferrite  material.  The  dispersion  relation  and  group  time 
delay  characteristics  of  the  waves  are  investigated. 

Magnetostatic  waves  are  slow  dispersive  magnetically  domi¬ 
nated  spin  waves  that  propagate  in  magnetically  biased  ferrite  slabs 
at  microwave  frequencies.  The  most  common  lowloss  ferrite  material 
used  for  MSW  propagation  is  epitaxial  yittrium  iron  garnet  (YIG). 

The  recent  interest  in  MSW  devices  at  microwave  frequencies 
has  accelerated  because  the  growth  of  uniform,  high-quality,  low- 
loss  epitaxial  YIG  films  with  a  large  aspect  ratio  has  been  improved 
and  the  design  and  fabrication  of  efficient  RF  to  MSW  transducers 
have  been  developed  and  realized.1  The  first  development  provides 
a  uniform  internal  dc  magnetic  bias  field  so  that  many  inhomogeneous 
transmission  problems  associated  with  nonuniform  internal  fields  are 
eliminated.  The  second  development  reduces  coupling  losses  from  RF 
waves  to  magnetostatic  waves  so  that  the  overall  insertion  loss 
can  be  reduced  greatly.  Another  motivation  for  exploring  the  use 
of  magnetostatic  waves  in  microwave  signal  processing  is  due  to 
the  fact  that  MSW  operation  is  possible  in  the  frequency  range  of 
1.0  to  20.0  GHz  with  wide  instantaneous  bandwidths  in  contrast  to 
conventional  devices  such  as  surface-acoustic-wave  (SAW)  devices 
which  usually  operate  on  the  IF  signal. 
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A  comparison  of  the  relative  merits  of  MSW  devices  and 
SAW  devices  is  given  by  Owens  et  al.2  and  Collins  et  al.3  and  is 
presented  in  Table  1.1.  From  this  table  it  can  be  seen  that  the 
tunable  properties,  the  lower  propagation  losses  at  microwave 
frequencies,  the  invariant  transducer  geometry,  and  the  adjustable 
delay  properties  of  magnetostatic-wave  devices  present  a  major 
advantage  in  terms  of  device  performance  over  other  existing  devices 
such  as  SAW  devices. 


Table  1.1 


Characteristic  Features  of  MSW  and  SAW  Devices 


2  ,3 


Property 


MSW 


SAW 


Maximum  delay 
per  cm 

Transducer  geometry 

Transducer  dimensions 


Hundreds  of  ns 


Simple /moderately 
complex 

Invariant  with 
frequency 


Hundreds  of  ps 

Complex 

Decrease  as  fre¬ 
quency  increases 


Table  1.2 


Prime  Applications  of  MSW  Devices 


Device 

Microwave  Function 

Nondispersive 
delay  line 

Phase  locking  of  pulsed  oscillator, 
signal  correlation,  communication  path 
length  equalizer,  rate  sensor 

Dispersive 
delay  line 

Group  delay  equalizer,  pulse  compression, 
compressive  receiver,  rate  sensor, 
frequency  synthesis 

Tapped  delay 
line 

ECM  deception,  PSK  matched  filter,  Fourier 
transformation 

Variable  delay 
line 

Target  simulation,  electronic  timing 

3andpass  filter 

ECM,  radar,  communication  satellite 
repeaters 

Tunable  resonator 

Narrow-band  frequency  filter  and 
oscillator  applications 

MSW  directional 
coupler 

Signal  routing,  switched  delay  line 

MSW  oscillator 


Stable  microwave  source 


microwave  communication  systems.  As  noted  in  this  table,  their 
prime  applications  are  in  the  areas  of  delay  lines,  filters, 
oscillators,  and  resonators.10  19  The  conventional  means  of  exciting 
magnetostatic  waves  are  by  metal  transducers  that  have  been  analyzed 
extensively  in  the  literature.20  35 

1.2  Previous  Investigations. 

Previous  analyses  of  various  types  of  geometries  for  MSW 
propagation  have  included  the  following: 

1.  The  first  structure  consisted  of  a  YIG  film  deposited  on  a 
substrate  (see  Fig.  1.1)  and  was  investigated  by  Damon  and  Eshback.36 
The  input  and  output  transducers  were  microstrip  lines.  They  deter¬ 
mined  the  different  propagating  modes  that  can  exist  in  this  structure 
A  summary  of  all  the  propagating  modes  for  different  magnetization 
directions  was  reported  by  Adams  et  al.9 

2.  A  multilayer  planar  structure  with  ground  planes  was  con¬ 
sidered  by  Tsai  et  al.37  and  others  ’  (see  Fig.  1.2).  In  this 
work,  wave  propagation  in  a  normally  magnetized  structure  of  infinite 
width  was  analyzed.  Daniel  et  al.40  reported  a  complete  summary  of 
the  wave  propagation  in  this  structure  for  principal  directions  of 
magnetization. 

3.  Young41  considered  wave  propagation  in  a  metallic  trough 
partially  filled  with  YIG  material  as  shown  in  Fig.  1.3.  He  derived 
the  dispersion  relation  for  two  important  propagating  modes  that 
can  exist  when  magnetization  is  parallel  to  the  slab  plane. 

U.  Finally,  a  long  rectangular  YIG  rod  with  metal  boundaries 
was  investigated  by  Auld  and  Mehta42  (see  Fig.  l.U).  They  studied 


FIG.  1.2  YIG  SLAB  BETWEEN  TWO  GROUND  PLANES. 
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two  possible  modes  of  propagation  in  this  structure,  using  a  mode 
analysis  technique,  and  derived  a  dispersion  relation  for  each  case. 

It  is  to  be  noted  that  in  all  the  reported  structures,  the 
ferrite  is  inside  an  unbounded  space  except  when  the  guide  is  com¬ 
pletely  filled.  The  case  of  a  ferrite  slab  enclosed  in  a  waveguide 
as  shown  in  Fig.  1.5  has  never  been  studied.  This  study  will  con¬ 
centrate  on  the  analysis  of  this  general  structure. 


1.3  Outline  of  the  Present  Study 

The  purpose  of  this  study  is  to  analyze  magnetostatic-wave 
propagation  in  YIG  slabs  placed  inside  a  rectangular  waveguide  as 
illustrated  in  Fig.  1.5.  The  general  formulation  of  the  equations 
for  magnetostatic-wave  propagation  and  some  known  results  for  cer¬ 
tain  special  structures  are  given  in  Chapter  II. 

In  Chapter  III,  the  case  of  a  ferrite  slab  in  contact 


with  both  sidewalls  of  the  waveguide  is  discussed.  This  part  can 


be  treated  as  a  boundary  value  problem  and  the  method  of  mode 
analysis  is  conveniently  applied.  In  this  method  a  certain  mathemat¬ 
ical  fcrni  for  the  potential  function  in  each  region  is  adopted  so 
that  it  satisfies  the  boundary  conditions.  Matching  the  potential 
fields  and  normal  field  components  at  the  slab  interfaces  gives  a 
system  of  linear  equations  which,  upon  mathematical  manipulation, 
yields  the  desired  dispersion  relation. 

When  the  width,  of  the  slab  is  less  than  the  v'dth  of  the 
guide,  that  is,  when  #  0  (Fir.  1 .  ‘  ) ,  the  mode  analysis  appears 
to  be  fruitless  and  ti.e  into. -ml  e  pint  h  r.  sec::..:  *.u  be  more  appropriate 


For  simplicity,  the  slab  is  assumed  '  t 


.In  it  ►  1)1  and 


s- 
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mathenatically  it  represents  the  source  of  magnetic  charges.  In 
this  way  the  integral  equation  is  tractable  and  yields  the  disper¬ 
sion  relations.  This  is  given  in  Chapter  IV.  Parallel  and  normal 
magnetization  are  considered  separately  in  the  wave  analysis  of 
Chapters  III  and  IV. 

In  Chapter  V,  numerical  methods  employed  in  obtaining 
numerical  data  and  dispersion  plots  are  presented.  The  results 
presented  in  this  chapter  extensively  describe  the  device  behavior 
and  performance  under  different  geometrical  configurations. 

Finally,  Chapter  VI  summarizes  this  study,  sets  out  con¬ 
clusions,  and  offers  suggestions  for  further  study.  Appendix  A 
contains  sample  computer  programs  developed  to  investigate  the 
structure  shown  in  Pig.  l.r  for  II  in  the  x-  and  z-directions. 


CHAPTER  II.  MAGNETOSTATIC-WAVF  PROPAGATION 


2.1  Introduction 

In  this  chapter,  the  mathematical  foundation  for  mapnetostatic- 
wave  propagation  in  unbounded  and  bounded  ferrite  media  is  introduced 
and  the  governing  equations  are  derived.  With  the  help  of  a  per¬ 
meability  tensor  derived  for  an  anisotropic  magnetic  medium,  it  is 
shown  that  plane  wave  propagation  of  electromagnetic  waves  in 
unbounded  ferrites  leads  to  three  regions  of  interest  in  the 
frequency  (w)-wave  number  (KJ  plane.  It  is  seen  that  magnetostatic- 
wave  propagation  is  possible  only  in  a  limited  range  of  wavelengths. 

In  this  range  of  wavelengths.  Maxwell's  equations  can  be  simplified 
which  makes  it  possible  to  ignore  the  electric  field  and  derive  the 
magnetic  field  directly  from  a  scalar  potential  field.  When  this 
approximation  is  used,  tic-  govern i n-'  partial  differential  equations 
inside  the  YT1  slat  ran  derived. 

Wave  propagation  in  a  lea  :e<:  rectangular  waveguide  under 
the  magnetostatic  at  prvxima4  :  •  >.  i.  .Jiscusr-ed  and  a  set  of  boundary 
conditions  which  mu  ft  he  satisfied  ai  the  metal  surfaces  and  slab 
interfaces  are  introduced,  finally,  sore  ki.'vn  structures  are 
analyzed  and  the  results  are  presented. 

2.2  Derivation  of  Effective  Permennil ity . 

Ferrites  are  ceramic-:  ike  magnet  i  <•  materials  where  high, 
resistivity  mane.;  <•  :dy  curren’  ._■%:>  remely  low  at 

microwave  frequencies.  The  *  cr,;-,rr  t  •  emeu  t  :  1  i  ty  of  a  ferrite  is 
derived  using  Hew*-!  ’  :•  ■  qua*  ;  •  f  t  i  where  t  iip  exchange- 


forces  ar.d  core  losses  are 


fif  =  ,  (2.1) 

where  M  =  the  total  magnetic  moment  i  er  unit  volume, 

It  =  the  total  internal  magnetic  field  and 

Y  =  the  gyromagnet ie  ruti o  which  is  equal  to  c/..e,  in  Cits  system 

Assuming  harmonic  time  dependence  Eq.  2.1  may  be  expressed  as  follows: 


O 


where  M  is 

n 

field  vector 


+  me' 


/—  ^  —  .not 

t.n  +  r.e 
o 


the  saturation  magnet. i nation  vector;  H  is 

o 

’  inside  the  slab  which.,  due  to  the  neglect 


)]  ,  (2.2) 

the  ac  magnet 
of  demagnet i- 


ic 


zafion  fields,  is  the  same  as  that 
are  the  RF  magnetic  moment  per  uni 
intensity  vectors,  respectively. 

Assuring  the  dc  fields  are 
neglecting  the  terms  of  order  high. 


‘‘  6  4  7  — 

outside  the  slab;  and  m  and  1 

volume  and  the  RF  magnetic  field 

,uch  larger  than  the  RF  fields  and 
•r  than  the  firs'.  ,  ho,.  2.2 


becomes : 


duor  =  u  y  (M  *  h.  +  m  *  H  ) 

c  c  o 


if  H  is  it:  the  z-di recti on  ar.d  h  is  nerren.iicu.lar  tr  H  ,  then  ho. 

i  ‘  go 

,'i  "  (2CiI'i  LC  OX’  1*0 S SCO  t  T;  "td’TTiS  C  T  i  "*  c C  V. '  f  g  t\y  '  1  "•  ■  vl'  ! 


Expressing  these  results  in  terms  of  a  permeability  tensor  yields 


and  the  free-soace  wave  number,  respectively ,  and  H  is  the  internal 

-  ex 

exchange  field.  Maxwell's  equation  in  conjunction  with  the  permeability 
tensor  from  Eqs.  2.1*  can  be  used  to  study  plane-wave  propagation  in 
unbounded  ferrite  media. 

2.3  Plane  Wave  Propagation  in  Unbounded  ferrites 

Magnetostatic  waves  have  wavelengths  much  greater  than  the 
lattice  spacing  and  therefore  it  is  appropriate  to  use  classical 
theory  rather  than  quantum  theory.  The  small-signal  theory  of  a 
lossless  ferrite  is  based  on  Maxwell’s  equations.  They  can  be 
written4'5’50  in  a  form  that  includes  the  permeability  tensor  of  the 
ferrite.  A  plane  wave  solution  of  the  form  ^  is  assumed. 

It  therefore  follows  that: 

?  *  h  =  jwse  , 

V  x  e  =  -  *i  :j  u  M  ■  h  ,  (2.?) 

o 

V  •  u  =  U 

and 

V  •  ( ,  .  •  h )  =  ii 

Solutions  to  Eq.  2.5  were  discussed  by  Seidel.51  It  is  con¬ 
venient  to  assume  that  the  magnetic  bias  field  is  in  the  z-direction 
and  the  propagation  vector  is  in  the  y-z  plane  at  an  angle  6  with 
the  7.  axis.  Including  tic  •  xchfinge  term  in  the  permeability  tensor, 
the  dispersion  relation  Is  found  to  be: 


-16- 


2  (u2  -  ii  -  K2  )  sin2  0  +  2?u  ±  !  ( 


H2  -  y  -  K2)  sir.14  6  +  UK2  cos2  e]2 


2! (u  -  1 j  sin2  e  +  T] 


,  (2.6) 


where  K2  =  u.'2ep  .  The  corresponding  RF  electric  and  magnetic  fields  are: 
o  o 


p(K  /K): 
o 


_ cos  e _ 

sin2  0  -  (K  /K)2 


K  tv  / 
or  o 


cos  e 


-  sir.  6 


-  ,1 K  •  r 


(K  /K ) ‘ 
o 


u(K  /K)2  -  1 
o 


-JK*r 


(2.8) 


p(K  /K)s 
c 


cos  0  sin  6 
*  sit,7  0  -  (K  /K)2 


Equation  2.6  defines  a  manifcl.a  of  dispersion 


tour:  dec 


by  curves  0  =  0  degrees  and  0  =  90  degrees  as  shown  in  'Fig.  2.1.  At 
each  value  of  0  the  medium  supports  two  types  of  waves,  an  ordinary 

wave  for  which  K/K  s  1  and  an  extraordinary  wave  for  which  K/K 

o  ‘  o 

becomes  very  large  at  certain  frequencies.  As  seen  in  Fig.  2.1, 
the  inclusion  of  exchange  forces  causes  the  lower  branches  of  the 

extraordinary  curves  to  bend  upward  for  large  values  of  K/K  . 

o 

The  lower  branches  in  Fig.  2.1  may  be  divided  conveniently 
into  three  regions. 

1.  K/K  ?  1.  Here  the  electric  and  magnetic  fields  of  Fqs. 

2.7  and  2.8  are  of  comparable  magnitude, 
region  of  "electromagnetic  prorogation." 


Hi  is  may  be  termed  the 


ORDINARY  WAVE 


2.  1  <<  K/K  <<  u)  /{•>,  i‘K‘).  i-rom  Eq.  2.7  it  is  seen  that 

o  o  ex  o  ^ 

the  electric  field  is  negligibly  small  in  this  region  compared  with 
the  magnetic  field  and  can  eventually  be  neglected.  Equation  2.8 
becomes : 


h 


-  J/K, 


-,iK-r 


(2.9) 


where  terms  of  order  (K  /K)^  have  been  neglected.  It  is  shown 

o 

easily  that  the  magnetic  field  in  Fq.  2.v  satisfies  the  magnetostatic 
equations  which  can  be  derived  from  Eq.  2.5  as  follows: 

V  x  h  =  0  (2.10a) 


and 


V  •  ( p  •  h )  =  0 


(2.10b) 


This  may,  therefore,  be  termed  the  region  of  "magnetic  propagation. 
It  corresponds  to  the  very  flat  portions  of  the  dispersion  curves 

where  the  frequency  is  m2  =  j  Ui  +  s in 2  0).  These  waves  are, 

o  o  M 

therefore,  characterized  by  a  low-phase  velocity  and  a  negligibly 
small  group  velocity. 

3.  a:  /(w  P2K2)  <<  K/K  .  Here  the  magnetic  field  is  again 

o  ex  o  o 

given  by  Eq.  2.9.  However,  the  exchange  term  now  has  a  significant 
effect  on  the  shape  of  the  dispersion  curves  which  bend  upward  into 
the  familiar  exchange  spin-wave  manifold: 


(2 


0.  f 
ex 


2K2)( 


;in‘- 


0 ) 


This  may  be  termed  the  region  of  "exchange  t  reparation." 


These 
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regions  are  of  course  not  shortly  lefiue.i,  :•  ince  the  propagation 
characteristics  vary  continue -ugly  w i -  h  K/K  . 

Phe  wave  p  ro:  agatio;.  hoc, use:  here  was  in  the  absence  of  any 
surfaces  an:  basically  ignore.:  t he  finite:  <iin.enr.i<  i.s  of  the  ferrite 
material  which  is  of  significant  value  to  engineering  problems.  Most 
of  the  structures  studied  so  far  in  the  literature  are  slabs  of 
finite  thickness  and  infinite  width.  In  the  next  section  +he 
governing  equations  for  magnetostatic-wave  propagation  are  derived 
and  the  boundary  conditions  that  must  be  satisfied  by  the  propagating 
waves  are  given.  Finite  dimensions  of  the  ferrite  material  and  the 
relative  orientation  of  the  magnetic  bias  field  with  respect  to  the 
propagation  vector  K  lea i  to  different  modes  of  propagation  which  are 
discussed  ir.  the  next  section. 


:;.l  Magnetostatic-V.’ave  S  rci  agat  i<~r. 

Wave  propagation  ir.  a  ferrite  medium,  of  infinite  extent  was 
considered  in  the  prey:  .us  sect  ion .  The  val id : 1 y  of  this  idealization 
as  an  approximation  tc  the  magt.*»t  ost  ati  c-vave  r  agation  in  a  finite 
medium  oi  a  rectangular  waveguide  is-  examine..:  next  ( pee  1.5). 

I'l.e  transverse  dimer. si-  ns  >  f  the  va.vgu  nr>-  sr.ai  i  e.-mi  a  red  to  the 


ileetromagnetic  wave  len-'th  am  Iht-  >■'.  eo*  r  rt.'i  iv 


■u'c  either 


cut-off  waves  or  leaky  wave r . 


2 . ‘i .  i  Govern  i  r.g  r-ount  . 


was  .-.hrwri  that 


Maxwell's  equations  it. 


s.tnt  !>•  a:  ;  ruximat  i  r-r.  reduce  to: 


'hut; 


i 


This  permits  the  definition  o1'  a  magnetic  scalar  potential  $  by- 
defining  h  =  V<f>.  The  magnetostatic  field  inside  the  ferrite  is 


governed  by  the  equation: 


pj  V  =  0 


.2.121 


When  H^c  is  either  in  the  x-  or  z-direction,  Eq.  2.12  takes  the 


following  forms: 


and 


>  +u(<f>  +  4>  )  -  0  ,  H  II  z 

zz  H  Yxx  Yyy  ’  dc  11 


t>  +  u  ( <{>  +  <(>)  =  0  ,  H  II  x 

xx  p  Yzz  Yyy  dc  11 


(2.13a) 

(2.13b) 


The  scalar  magnetic  potential  satisfies  the  following  equation: 


or 


outside  the  ferrite. 


V^<f>  =  0 


f>  +  <{,  +4,  =0 

xx  yy  rzz 


(2.13c) 


jwt 


In  all  succeeding  analyses >  the  time  harmonic  variation  e 

is  assumed  and  is  omitted  for  simplicity.  The  wave  propagation  is 

—  i  Kv 

in  the  y-direetion  and  varies  as  e  .  Thus  Eqs .  2.13  reduce  to 


XX 


rXX 


in  the  YIG 

slab: 

+  <{> 

zz 

mk2  4> 

’  Hdc 

+  u<|>  = 

zz 

V>K  2  (j> 

•  Hdc 

and  in  ai 


<p  +  f  -  K‘ 
xx  fzz 


(2.ll+a) 

(2.ll+h) 

(2.lkc) 


2.b.2  Boundary  Conditions.  The  YIG  slab  is  enclosed  in  a 
metallic  waveguide  and  the  following  boundary  conditions  should  be 
satisfied  at  all  metallic  guidewalls  (see  Fig.  1.5): 


X 


(2.15a) 


It  t  M 


=  0  i 4 


(2.15b) 


As  seer:  in  iir.  1.5,  *  he  o:  :i-..iti  :r;3  t  hut  must  be  satisfied  at 
the  four  slab-air  interfaces  are  as  fed]  ws : 


and 


3.  <t>  continuous  at  z  7  , 

U.  $  continuous  at  x  =  x  ,  (a  _  v  ) 

o 

5.  b  continuous  at  z  =  z.  ,  z, 

Z  12 

6.  bv  continuous  at  x  =  x  ,  (a  -  x^) 


(2.16a) 


(2.16b) 


(2. l6c ) 


(2.1bd) 


In  general ,  this  pro bier:  is  lest  solved  by  the  integral  equation 
method  where  the  slab  can  be  considered  to  be  the  source  of  magnetic 
charges.  This  formulation,  though  complex,  is  very  successful  when 
■  t.e  slab  width  is  less  than  the  guide  width  and  provides  accurate 
answers  to  the  dist-eruir.r.  cal cu'.at ionu  required  for  plotting  the 
curves .  This  general  case  is  discussed  and  analyzed  in  Chapter  IV. 
d:  wover,  in  the  spec:.':]  ..-are  when  the  Yid  slab  extends  to  both  guide 
sidewalls,  the  prcbl  om  ear  t  re-,  t  ed  ns  a  boundary  value  problem 
and  a  simpler  method,  i.e.,  new-  analysis,  w.uli  b*»  very  effective 
1:;  this  case.  This  case  is  considered  ir.  Chapter  111. 

2.b.3  Modes,  of  Propagation.  Magnetostatic-wave  •  reparation 
in  thin  fi  lms  has  been  invert,  igul  e  i  ox*. « ••.rive!;.  and  *-hroo  major 


5  t.  -  f  o 


i  ropapati nr  modes  have  been  usv*d  :  n  :ev:ce  art  1 1  cat  i  c  no  . 


These  three  modes 


bias  magnetic  field,  propagation  direct i  u. ,  and  the  fvrr i t r  s  an. 
These  modes  are  disnerrive  and  '•'.aruc*  crimed  t  y  a  1  :r.:  *  •  :  :  r  :  urnt  ior. 
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passband  and  magnetic  bias  tunability.  The  following  is  a  brief 
discussion  of  these  modes. 

1.  Tiie  magnetostatic-surface-wave  (MSSW)  mode  in  which  the 
bias  field  H  is  perpendicular  to  the  direction  of  wave  propagation 
and  both  are  in  the  plane  of  the  film.  This  mode  has  highly 
anisotropic  propagation  in  the  film  plane  and  the  mode  energy  is 
confined  to  the  "top"  surface  for  forward  propagation  and  to  the 
"bottom"  surface  for  reverse  propagation. 

2.  A  second  mode  is  the  magnetostatic-forward-volume  wave 
(MSFVW)  in  which  the  bias  field  H  is  perpendicular  to  the  film 
plane  and  is  characterized  by  approximately  isotropic  propagation 
in  the  plane  of  the  slab.  The  lowest-order  mode  is  usually  most 
easily  excited. 

3.  The  third  mode,  the  magnetostatic-backward-volume  wave 
(MSB VW)  exists  when  the  direction  of  the  bias  field  and  propagation 
are  the  same  and  in  the  film  plane.  This  mode  has  opposite  phase 
and  group  velocity  directions,  is  highly  anisotropic,  and  is 
multimoded  as  in  the  MSFVW  case. 

Figure  2.2  shows  the  dispersion  characteristics  of  the  three 
modes  of  operation  propagating  only  in  certain  frequency  ranges.9 
In  this  figure  the  dispersion  characteristic  of  a  surface  wave  when 
one  of  the  surfaces  is  metallized  is  also  shown .  As  can  be  seen, 
t  he  r.rcpagati  n  bandwidth  increases  and  device  performance  is  now 
iif foment  . 

The  idea  f  ;  iaeing  ground  $  lane;-  in  the  vicinity  of  a  YIG 
slab  , ana  ♦  bus  ■•hanging  i  *  .•  dirt  <•!>'•  :v*-  characteristics  and  its  group 

7  - o 

del  ay  was  i  nv  gt  i  ,-ti‘  •  d  .  i  r:  t  ho  ft!  1 w  5  nr  section,  the 


MSBVW 


progressive  structural  development  3f  the  device  is  cited  and  briefly 
analyzed.  In  most  t>f  there  investigated  structures,  the  YIG  slab  is 
assumed  to  have  finite  thickness  and  infinite  width.  The  latter 
assumption  is  not  realistic  in  -a.  engineer: ng  sense  and  needs  to  be 
modified.  In  Chapters  III  and  IV,  finite-width  slabs  are  considered 
under  general  boundary  conditions. 

2.U.U  Results  of  Previous  Investigations.  A  literature 
review  of  propagating  waves  in  several  structures  is  presented  next. 

1.  Unshielded  slab  in  free  space.9  This  configuration  is  shown 
in  Fig.  2.3  where  three  principal  directions  of  magnetic  bias  fields  are 
shown.  The  three  different  cases  are  discussed  next. 

a.  H  in  the  x-direction.  The  propagating  waves  in  this 
case  are  called  magnetostatic-surface  waves  (MSSW)  and  the  dispersion 
relation  is  given  by:61 


2Kt 

e 


K/2)2 


UJ  2)]2 


(2.17) 


where  t  =  the  slab  thickness. 

b.  H  in  the  y-direction.  The  prODagating  waves  in  this 
dc 

case  are  termed  magnetostatic-backward- volume  waves  (MSBVW)  and  the 
d i spers i on  relation  is  found  to  be ; “fit 


2  cot  aKt  =  a  -  a~ 1 


(2.18) 


where 


«x 


c.  H  in  the  z— direct?  on.  These  waves  are  called 
dc 

magnetostatic- forward- volume  waver,  { M2  FYW )  and  the  dispersion 


relation  is  given  by: 


f,  r, ,  6  r, 


tan  (Kt/2u) 


(2.19) 


where 


u  ('i)  +  -  ,*'2 

o  o  M 


The  three  waves  are  for  a  slab  extending;  in  width  and  length  to 
infinity  but  with  a  finite  thickness.  Equations  2.18  and  2.19 


provide  many  modes  due  to  the  sinusoidal  nature  of  the  dispersion 
relations . 

3  7  —  3  9 

2.  YIG  slab  between  ground  planes.  Two  cases  of  interest 

are  dicussed  (see  Fig.  1.2): 

a.  Farallel  magnetization.  The  dispersion  relation  is 
given  bv:40 


2Kt  (u  -  Ki  -  tanh  Ks^Xp  +  Kj  -  tank  Ks ?) 

(p  +  Ki  +  tanh  Ks  )(p  -  +  tanh  Ks^) 


!  2 . 20 ) 


b.  Normal  magnetization.  The  dispersion  relation  for  this 


care  is  found  tc  be; 


tar.  (oKt)  =  — 


tanh  K’*  tanh  K's  ’ 


(2.21 ' 


wh.ere  a2  =  -  I1  '' 

3.  Partially  oh  i  •  •  i  <led  slab  w  it.h  farallel  magnetization.  The 
geometry  for  this  ease  is  sh-  vn  in  1  2.4.  Two  modes  of  propagation 

art-  possible  in  4  hi.  ms-  ’..'.s:;  are  discussc-i  below. 

■i.  Pl.t  first  aris's  ’.lie  configuration  of  Fig.  2.4a 


and  is  c.r: 


r<  f  1  f,  j  4  !;<•  ferri  *  t— ::.ofa  1  P:M 


The  surface  wave  clings 


to  the  shielded  r.ur  fa  The  rs : ersi  or  rei'i+Pn  is  given  by 


r.  = 


Yn/2Kl){tl 


"W 


1  ( my  /y 1 

n  r  i  ■ 


:ot'r 


Y  a 
n 


-  1} 


(2.22) 


where 

Y  =  {[  (Wfc)2  +  uK?  Vm  }  * 


n  =  0  ,1 


and 


=  f  ( m  /b ) 2  +  K 


2  ' 


b.  The  second  mode  arises  in  the  configuration  of  Fig.  2.«b 
and  is  referred  to  as  the  ferrite-air  (FA)  ir.cde.  Here  the  surface 
wave  fields  cling  to  the  unshielded  surface.  The  dispersion  relation 
is  given  by41 

K  =  (y^KjHU  +  ^(py^/yM''  +  “(uY„/y',)  =  oth  Yr,a’r  +  1}  , 

r.  =  f,l,f  (2.2V 

Vlier.  the  direction  f  ;  ropagatior.  or  1  hr  magnetic  fiel  d  tires*  it:;  ar< 
reversed,  the  properties  of  the  FA  and  -V.  modes  art.  interchanged. 
These  are  shown  in  Fig.  2.H. 

A.  Completely  filled  guide.  The  struct ur*-  f.  r  this  .  •«!.*.•• 

is  shown  in  Fig.  ''.fa.  It  is  four.:  that  t  he  pro-}  ?igat  :  ::.*  ,t.  :-uvr  are 
either  surface  cr  volume  waves  that  :  re;  agate  in  :  if  fere:. t  f  re  luency 
ranges . 

fi •  V'y.iu:.’’"  rr/‘!:  i »  ’  ko  ••!.?  rr-err  "i  'n  ?v.I  Hf  i  r-u  *'  r  4:h:.' 

is  riven  by  ( n  ,:r  ^  1  ■ 


t.r  /t. ) '  /u 


(2.2V 


where  ,n  are  nonr.vr 


(.*  t  r  e  v . : 


H=250  Oe 
o-  0.0856  cm 
b  =  0.785  cm 


FREQUENCY 


b. 


given  by 


Surf ace  waves . 


She  dispersion  relation  for  this  case  is 


K2 


( nn /b ) 2 
(K2/u)  -  u 


U  >  0 


(2.25) 


These  two  cases  are  slotted  in  Fig.  2.6b  for  H,  =  500  Oe  and 

dc 

a  =  b.  In  this  particular  case,  the  dispersion  relation  can  be 
shown  to  be  independent  of  dimension  a.42  As  noticed  here  in  almost 
all  the  structures,  magnetostatic-wave  propagation  was  studied  in 
unbounded  media,  thus  giving  only  an  approximately  answer  to  the 
problem  of  the  wave  propagation  in  finite-width  YIG  slabs  in  the 
presence  of  boundary  surfaces  made  of  metal.  Under  these  general 
conditions  the  solutions  offered  here  by  the  widely  known  boundary- 
value-problem  technique  would  no  longer  hold  and  a  more  powerful 
technique,  i.e.,  the  integral  equation  method,  must  be  employed 


in  order  to  obtain  reasonable  answers. 


CHAPTER  III.  MODE  ANALYSIS 


3.1  Introduction 


In  this  chapter  magnetostatic-wave  propagation  inside  a 
waveguide  partially  filled  with  a  YIG  slab  (see  Fig.  3.1)  is  investi¬ 
gated  by  the  mode  analysis  technique.  The  slab  extends  to  both 


guide  sidewalls  and  therefore  the  problem  can  be  treated  as  a 
boundary  value  problem. 

The  differential  equations  to  be  solved  are: 

V  •  ( yV<j) )  =  0 

in  the  YIG  and  the  Laplace  equation  V 2 <J>  =  0  in  air.  The  potential 
field  in  the  air  or  YIG  regions  can  be  written  in  terms  of  a  Fourier 
series  expansion.  The  Fourier  series  expansion  in  each  region  should 
satisfy  the  boundary  conditions  at  all  metal  surfaces  as  given  in 
Section  2.1*.  2. 

Matching  the  potential  functions  at  the  slab  interfaces 
provides  a  system  of  linear  equations.  Rewriting  the  system  of 
linear  equations  in  matrix  form  and  requiring  a  nontrivial  solution 
to  the  problem  results  in  setting  the  determinant  of  the  coefficient 
matrix  to  zero.  Solving  this  equation  for  roots  yields  the  disper¬ 
sion  characteristics. 


3.2  Potential  Fields  in  the  Air  Region 


In  Fig.  3.1,  the  differential,  equation  that  should  be  solved 
in  the  air  region  is  the  Laplace  equation: 


Lx  +  <f>  ~  k2<!>  • 

XX  Z 


(3.1) 
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Solutions  to  Eq.  3.1  in  the  air  regions  (I  and  III)  can  be  written 


as  a  combination  of  even  and  odd  sinusoidal  harmonics  in  the  x- 
direction  and  in  terms  of  hyperbolic  functions  in  the  z-direction 
as  follows: 
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where 


r 


Equations  3-2  should  new  saiirfy  (>..  ui.  j*j ry  c  river,  in 

Section  2.U.2.  Fron;  Jojuiit  i n\  ,  l  =  f  ^/:.ix )  =  C  at  x  =  C,  a, 

Tliis  p;ives 


and 


r 

?■  n 


From  Condition  2,  b  =  ) 


This  tf’ives 
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Theref ore  ,  Eqs .  3.2s impl  i fy  *  c 


co 


on 


n=o 


(3.3) 


(3.U) 


where  A  and  D  are  unknown  constants  which  will  be  determined  when 
n  n 

the  potential  functions  and  norraal  magnetic  fields  are  matched  at 
the  slab  interfaces. 


3-3  Potential  r'ields  in  the  YIG  Region 

The  direction  of  the  bias  magnetic  field  in  the  cross-section 
plane  can  be  chosen  to  be  either  parallel  or  normal  to  the  YIG  slab. 
In  each  case  the  partial  differential  equation  for  the  potential 
fields  in  the  YIG  given  by  V *  ( p V <+> )  =  0  is  dependent  on  the  direction 
of  magnetization.  Therefore  each  case  is  analyzed  separately. 


1.  Parallel  magnetization.  The  relative  permeability  tensor 

when  H .  is  in  the  x -direction  (see  tig.  3.2)  from  ha.  2.*‘b  is  written 
n  c 


(3.6) 


'.>7. 


The  following  form  Is  a  s  dir  *  _•  •  nr.  i  satisfies  boundary 

ditioii  1  (Section  ;  ’ .  t . ;•  ir.  Y  : v  (II): 


oo 


r>=o 


where  B  and  C  are  arbitrary  constants  nnd 
n  n 


Y 


n 


(3.8) 


2.  Normal  magnetization.  The  do  magnetic  field  is  in  the 
z-iirect ion  (see  Kir.  3.3)  and  the  relative  permeability  tensor 
from  r’q.  2.Ua  is  given  by 


(3.9) 


,*rr  .at  ion : 


(3.1.  c 


lYi 


com 
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Equation  3 .  lU  represents  the  peter.'  1 


slab  in  terms  of  the  unknown  cons!  . 


potential  functions  at  the  slat  :i. 


relation  and  allows  the  1  etc- it. ’ 


3.U  Dispersion  Relation 

To  obtain  the  dispersion  r-1  a';  n, 

3  and  5  from  Section  2.b.2  must  be  a;:]  ; : 


magnetization  case  as  follows: 


1.  Parallel  magnetization.  Condition  3  r*.  :ulrv; 
of  the  potential  functions  3.3  and  3.1*  (in  air)  an  i  3." 


z  =  z  ,  z  .  This  gives 
12 


V  A  cos  —  x  cosh  v'(b  -  z  ) 
L  n  a  n  2 


-.IKy  _ 


y  cos  —  x(B  cosh  y  z„ 
L  an  'n  2 


+  C  si n'n  Y  2  ) 
n  'n  2 


,-JKy 


(3.15a) 


y  D  cos  —  x  cosh,  y  2  e  "  ' 
“  n  a  n  l 

n=o 


y  cos  —  X(B  cosh  y  2 

L  r>  1  n  1 


+  C  sinh  y  7  )  e 
n  '  n  l 


(3.15b) 


Condition  5  requires  continuity  of  b  at  the  interfaces 
z  =  z  and  z  which,  upon  equating  (from  Eqs.  3.3,  3.**  and  3.7),  gives 


y  -Ay*  cos  —  x  sir.:.  v'(b  - 
u  n  n  n 


+  uv  cc-sh  r  r? J  + 


)  cos  “  x ,  h  (-  K  K  cosh.  Y 
"  a  n  l  n 

i.= 


1  sinh  +  p-f,  •'  r-h  7 r, 2 2 )  1  (3.15c) 
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The  existence  of  nont rival  solutions  for  A  ,  B  ,  C  and  D  requires 

n  n  n  n  H 

the  determinant  of  the  coefficient  matrix  in  Eq.  3.16  to  be 
zero.  This  yields  the  following  dispersion  relation: 


tanh  y  t  I  KTK? 

rj  i 


y ‘  -  '  1  tanh  y's  tanh  y's  -  K  Ky  (tanh  y's 

'n  'r.  'n  2  'n  i  l  rn  'n  2 


-  tanh  y^si  )  ]  +  yy^y^tanh  y^s?  +  tanh  y's  )  =  0  ,  (3. IT) 


where  t  is  the  thickness  of  „he  TIG  slab  and  s  and  s  are  the  thick- 

l  ? 

nesses  of  air  regions  Ill  i. 

An  alternate  form  of  Eq.  3-17  is  given  by 
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(3.18) 

Equations  3. IT  and  3.18  are  exact  dispersion  equations  which  determine 


the  frequency  range  of  propagation  of  each  mode  independent  of  the 
existence  of  tile  other  modes. 

These  equations  reduce  tc  simpler  forms  when  some  or  all  of 
the  metallic  boundary  surfaces  are  driven  to  infinity.  All  these 
cases  are  studied  under  the  heading  of  ierer.erate  cases  in  a  latc-r 


section . 

2.  Ilormal  magnetization.  Matching  the  potential  functions  an: 
normal  magnetic  fields  in  the  air  '  Eqc.  3.3  and  3.M  and  in  TIG 
{ Eq .  2.lU)  at  the  interfaces  :•  =  r.  uni  r.  on  the  basis  of  ^  n  iiticn 
3  and  5  giver,  in  et i--n  2J..P  ••  i vvr 
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As  can  be  seen  in  Eqs.  3.19,  all  inodes  except  the  zeroth  order  are 

coupled  and  one  cannot  exist  without  the  others.  This  phenomenon 

leads  to  mode  coupling  whicn  complicates  the  elimination  of  the  unknown 

constants  A  ,  B  ,  C  and  D  , in  order  to  find  the  dispersion  relation. 

’  n  n  n  n 

However,  the  following  procedure-  is  adopted  to  obtain  the  dispersion 
relation:  Multiplying  both  sides  of  Eqs.  3.19  by  cos  (mn/a)x, 

the  interface  from  f.  to  a  and  noting  that 


integrating  along 


r 


Upon  elimination  of  and  from  Eqs.  3.20  the  following  system  of 
linear  equations  is  obtained: 


BR  +CS  +P  V  (BR  +CS) 
m  mm  m  mm  ^  n  mn  n  mn 
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=  0  ,  (3.21b) 


where 


R  =  8  (y'  tanh  y's  cos  y  z  -  y  sin  y  z  )  , 
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P  =  -  1+K^Ka/im 2 


Equations  3.21  provide  an  infinite  number  of  linear  equations  in 


B  and  C  which  can  be  expressed  in  the  following  form: 
n  n 


a.  Filled  guide.  This  ease  corresponds  +o  s  =  0  and 
s2  =  0  and  was  presented  in  Case  4  of  Section  2.U.U.  From  Eq.  3.17 
the  following  is  obtained: 


(tanh  y  b)(K2K2  +  p2y2)  =  0 

n  l  n 


(3.23) 


Setting  the  term  in  the  first  parentheses  to  zero  in  Eq.  3.23 


gives 


tanh  (y  b)  =  0 


(3-2U) 


If  jj  <  0,  the  solution  to  Eq.  3. '2b  is  given  by 


.  nu 

Yn  J  b 


(3.25) 


Combining  Eq.  3.25  with  Eq.  3.8  yields  Eq.  2.2b  precisely. 

Setting  the  terms  in  the  second  parentheses  to  zero  gives 


Y2  =  -  K2K2/y  . 

'n  l 


(3.26) 


Combining  Eq.  3.26  with  Eq.  3.8  yields  Eq.  2.25. 

b.  Infinite-width  YIG  between  ground  planes.  If  a 

the  problem  reduces  to  that  of  Case  2  in  Section  2 . U . U .  Equation 
3.17  yields  the  same  equation  as  Eq.  2.20  upon  proper  substitution. 

c.  Inf inite- width  slab  in  free  space.  Letting  a  ■+  oo, 

s  ->■  oo  s  -*  «  reduces  the  problem  to  that  of  Case  1  in  Section 
12 

2.U.U  and  Eq.  3.17  in  this  case  corresponds  to  Eq.  2.17. 

2.  Normal  magnetization.  Here  one  simple  case  of  interest  is 
derived  and  the  other  cases  simply  follow  the  same  procedure.  Through 
analysis,  it  was  found  that  the  zeroth-order  mode  is  uncoupled  and 


corresponds  to  the  case  when  a  ->  ® ,  i.e.  ,  when  the  YIG  slab  is 


shown  in  Fig. 


-Un¬ 


placed  between  two  ground  planes.  1'his  structure  as  £ 

3.U  has  been  investigated  extensively.  7,40 

i  V>  derive  the  d  i  r.pe.T  i  on  relation,  tine  ■  let,  erir;  i  riant  of  the 

zeroth-order  matrix  (M  )  from  Eq .  <.22  must  be  set  to  zero  as  follows 
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Upon  substitution  and  simplification,  Eq.  3.27  becomes 
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where 


bunt  ion 


3.1*  YIG  SLAB  BETWEEN  TWO  GROUND  PLANES  (NORMAL  FIELD). 
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3.  By  means  of  the  "rcaf.net ie  sources"  and  appropriate  Green's 
function  an  integral  expression  for  the  potential  function  everywhere 
inside  the  waveguide  (including  the  ferrite)  can  be  written  as  follows: 


( x  ,y  ) 


i;  Pv(x' ,y ' ,z' )Gi(x,y,z)  dv' 
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YIG 

surfaces 


where  ^(x^jz)  the  potential  function,  is  the  Green's  function, 
(x,y,z)  is  the  observation  point  and  (x',y',z’)  is  the  source  point. 

1*.  Assuming  wave  propagation  in  the  y-direction,  the  y-variation 
of  all  functions  involved  in  this  study  is  therefore  of  the  form 
‘  .  In  this  manner  the  following  can  be  written: 
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Fir.-  4.1;.  The  l  our.' j  ary  emir  ion  cn  G(x,z)  is  that  its  normal 
derivatives  should  vanish  <  t.  the  four  sidewalls. 

Following  the  conventional  approach,67  G(x,z)  can  be  expressed 
as  follows: 

G  ( x ,  z  )  =  G  ( x  ,z  )u( z 1  -  z)  +  G  ^ ( x  ,  z ) U ( z  -  z')  ,  (4.4) 

where  G  (x,z)  and  G?(x,z)  satisfy  the  homogeneous  partial  differen¬ 
tial  equation: 


V  Z 

T — 5-]g(x,z)  -  K2G(x,z)  =  0  . 


(4.5) 


The  boundarv  conditions  for  G  and  G  are: 
”  1  2 


-  0  at  x  =  0  ,  x  =  a  , 


at  z  =  0 


=  0  at  z  =  b 


Formal  solutions  for  G  (x,z)  and  G  (x,z)  can  now  be  written  as 
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x  z )  =  y  cos  —  cosh  '  (b  -  ?.)  z  >  z’  ,  (4.6b) 
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A  cosh  v’z*  =  b  cosh  y'(b  -  z') 
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(M) 


2.  To  obtain  the  second  source  condition,  Eq.  U.3  is  integrated 
from  z'  -  e  to  z 1  +  e  where  c  is  a  very  small  number.  This  gives 


■<r,  i  id 


2 

a(  1  +  6 


nn  . 
cos  —  x' 
a 


(U.8) 


Differentiation  of  Eq.  h.C  and  substitution  in  Eq.  U.8  gives 


B  sinh  v '  (b  -  z')  -  A  sinh  y'?-'  =  - r 

n  n  ri  'n  y  '  a(l  +  6  ) 


From  Eqs.  07  and  U.9,  A  and  P  are  obtained: 
^  n  n 


nn  . 

COS  -  X  ' 

a 


(h.9) 


(l  *  So „> 


cos  —  x'  cosh  v'(b 
a  n 

sinh  y ' b 
n 


;U.10a) 


Y'a(l  +  6  ) 

n  oi  i 


cos  — 1  x '  cosh  Y 1 ; 

_ a _ 

sinh  v'b 


(U.lOb) 


Using  Eqs.  h.C  and  H.10  yields  the  Treen'r:  function  0(x,z)  from  Eq. 
l.H  in  tiie  final  form: 


_ ,  .  .  ,  r  7  nn  ,  nn 

0  ( :< ,  x '  ,z,z'  =  )  -  \ - — - 77  cos  —  x'  cos  —  x 

’  ’  ’  L  v'a(l  +  <*.  '  f.ji.i.  /  b  a  a 

n=c  ’n  on  'n 


CO. oil  y’(b  -  cnrh  y'z  U  ( y  1  - 

1  n  ri 


Y 1 a( 1  +  <5  )  sinh  y  fb 

n  on  n 


nn  ,  rni  ,  .  ,  .  , 1 

cos  —  x  cos  —  x  cosh  y  z  cosri  y  { b 
a  a  1  n  n 


z )  V  (  z  -  z 1  )  .  ( . 1 1 ) 


1 . 3  Equivalent  Sources 


In  order  to  be  able  to  solve  the  integral  equation  in  terms 
of  the  potential  inside  the  Y  Ki  slab,  the  magnetic  charges  must  be 
determined.  Mathematically,  the  source  of  magnetic  charges  can  be 
considered  to  be  the  YIG  slab.  The  mathematical  1'orn:  of  the  magnetic 
charges  in  terms  of  ip  strongly  depends  on  the  permeability  tensor  which 
itself  is  a  function  of  magnetic  l>inc  field  direction.  In  this  work 
parallel  and  normal  magnetization  are  studied  and  therefore  the 
magnetic  charges  are  derived  separately  for  each  case. 

To  derive  the  magnetic  charges,  the  small-signal  magnetiza¬ 
tion  vector  (m)  must  be  determined.  When  the  magnetostatic  approxi¬ 
mation  discussed  in  Section  2J1.I  is  used,  the  magnetic  field  intensity 
(h)  in  the  slab  can  be  derived  as  the  gradient  of  a  scalar  potential 
function  <j>(x,z ) .  Finding  k  in  terms  o*’  <j>(x,s)  leads  to  determination 

f  rti  and  eventual  ly  to  p  as  fol  levs': 

rn 

:•  ---  c  +  p  ,  (t.h.:) 

It. 

w:.  •  re 

Pv  =  -  7*m  +  Xrx.  f  ,  ( a .  1 3a ) 

(  =  ir.  •  n  ( .  1  M 1 

••it:  i  i,  is  a  unit  vector  normal  to  the  slab  surface.  The  rmal  1  -r i rrnl 
magnet  i  zatior.  (m)  is  rivet:  by 
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and  is  the  permeability  tensor  derived  in  Section  2.2.  Substitution 
for  b  and  h  in  Eq.  4.l4  gives 


m  =  (pr  -  I )  *  ( V <f>  -  jK<j>y )  ,  (4.16) 

where  I  is  the  identity  tensor,  i.e.,  all  elements  are  zero  except 
the  diagonal  elements  which  are  1. 

Thus  with  the  help  of  Eq.  4.l6,  surface  and  volume  charges  from 
Eqs.  I*.  13  can  be  written  as: 

Pv  =  -  V-[(ur  -  l)*(V(j>  -  jK^y )  ]  +  jKmy  (4.17a) 

and 

Ps  =  f(br  -  T) - (V*  -  jK*y)].n  .  (4.17b) 


From  Eqs.  U . 17  it  can  be  seen  clearly  that  the  magnetic  volume  and 
surface  charges  depend  upon  the  relative  permeability  tensor  (ur).  In 
the  following  sections,  and  pg  are  derived  explicitly  in  terms  of  the 
scalar  potential  in  the  slab  (<+>)  for  the  two  cases  of  parallel  and 
normal  magnetization. 

4.3.1  Parallel  Magnetization  Charges.  The  relative  permeability 
tensor  from  Section  2.2  for  H  in  the  x-direction  is  given  by  (see 
Fig.  4.2): 


1  0  0 

0  V  JK1 

0  -  ,iKi  y 


(4.18) 


Substitutions  from  Eq.  4.18  into  Eqs.  4.17  yields  p  and  p  as 

v  s 


follows : 
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-  V  •< 


(„  -  i)  *  KiK* 


X  + 


JKX  ||  -  JK(W  -  1)1 


+  KjK  -  K2(u  -  1)$  (4.22a) 


3x 


and 


+  < 


— 

_  — 

(„  -  1M4 

X  + 

'  jKl  ft  ■  jK(y  -  1}* 

_ 

y 

( 4 . 22b ) 


where  the  upper  and  lower  signs  are  for  charges  at  x  =  x  ,  a 

o 


respectively.  Upon  performing  the  mathematical  operations  in 
above  and  noting  that  (from  Section  2.4.1): 


o 

the 


324> 

7T  +  u 


3z" 


ax" 


-  K2  4> 


0 


Pv  and  pg  are  finally  given  by 


P  (x,z)  = 

v 


V  -  .1  324> 


3z" 


(x,z) 


(4.23) 


and 


P  (x,z)  = 


a^(xQ,z) 

)  =  -  (y  -  1)  - r -  -  K,K4>(x  ,z)  at  x  =  x 

s  ^  ox  1  o  o 


\p  =  (p  -  i: 


3$(a  -  x  ,z) 


3x 


+  KxK^(a  -  xq,z) 

at  x  =  a  -  x  .  (4.24) 


From  Eq.  4.24  it  can  be  seen  that  the  surface  charges  at  z  =  r  or 


are  absent  and  the  only  existing  surface  charges  are  at  x  =  x^  and 


x  =  a  -  x  .  It  will  be  seen  that  this  difference  in  the  charge 
o 


arrangement  and  mathematical  form  for  the  two  cases  will  lead  to 


different  formulations  in  each  case. 


As  discussed  in  Section  k 


the  integral  equation  (Eq.  4.2) 


can  be  obtained  in  terms  of  {  .x  ,z)  i cr.iy  the  reints  insi  de  the 
ferrite  slab  are  considered.  Since  c  and  ,  are  different  for  the 
two  cases  of  parallel  arid  normal  magnetization,  each  case  is 
considered  separately. 

4.4.1  Parallel  Magnetization  rormulat  ion.  When  the  appropriate 

expression  for  and  p?  from  Fqs.  it. 19  and  4.2 C  is  used  and  when 

points  (x,z)  are  considered  such  that  {(x,z),  xq  =  x  =  a  -  x^ , 

<  <  .  | 

z  =  z  =  z  ; ,  Eq.  4.2  can  be  written  as: 


<j>  ( x ,  z  )  = 


a-x  z 
r  O  r  2 


O  1 


[K2(y  -  l)|(x',z')  -  (p  -  !)<}>  (x',z')] 


ZZ 


a-x 


»G(x,x',z,z')  dx'  dz'  + 


a-x 


•G1(x,x',z,z1)dx'  + 


f  K^K^(  x  '  ,z^  )  -  (p  -  l)^(x»  z  )] 
KiK$(x’,z2)  +  (p  -  l)^(x',z2)] 


Gz (x,x' ,z,z2 )  dx'  .  ( 4.25) 


Trie  integro -differential  equation  (Eq.  4.25)  is  two-dimensional  in 
(x,z)  and  is  difficult  to  analyze.  Thus  the  ferrite  slab  is  assumed 
to  be  very  thin  such  that  $(x,z)  may  be  assumed  to  vary  linearly  in 
z.  In  this  approximation  two  functions  of  one  variable  may  be  used 
to  approximate  <j>(x,z)  in  the  following  manner.  If 

f  ( x )  =  a  ( x ,  z  ) 

1  i 

and 

f  (x)  =  <j>(x,z  ) 

2  2 

then  by  linear  approximation  the  following  is  obtains: : 
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4>(  x  ,z )  s  f'1(x)  + 


f  (x)  -  f  (x) 

2  1  ,  s  <  < 

-  (2  -  2  )  2  —  Z  —  Z 

2-2  11  2 
2  1 


f2(x)(z  -  Z1 )  +  fi(x)(z2 

_____ 


(4.2 6) 


which  gives 


f  (x)  -  f  (x) 

f>  U,z)  £  - 

2  Z  -  Z 

2  1 


(it. 27) 


.(x,z)  s  0 


(it. 28) 


Substitution  of  Eqs.  4.26  through  4.28  in  the  integral  equation  (Eq. 
4.25)  yields 


a-x  z 
r  0,2 


e  U  e  * 

4>(x,z)  =  K 2  (  u  -  1 )  - 

J  x  7. 


1  x  '  z 
o  1 


f2(x’ )(z'  -  zx)  +  f1(x')(z2  -  z') 
z  -  z 

2  X 


G  ( x  ,x'  ,z  ,z'  )  dz  '  dx '  - 


fj(x')(z'  -  z^)  +  f 1 (x  *  )  ( z  2  -  z' 
( z  -  z  ) 

2  X 


+  ( 11  -  1 ) 


fz(x')  -  fi (x ’ ) 

(z  -  z  ) 

2  1 


G  (x,x'  ,Z,Z1 


r  |-f?(x')(z'  -  Z1  )  +  f](x')(z2  -  z1 

L  v:i 


+  ( u  ~  3 


r  ( X ’  )  -  r  (x'Kl 


(z2  -  Zl} 


G ,(x,x' ,z ,z2 )  dx'  . 


(4.29 ) 


Evaluating  the  integral  equation  given  by  Eq.  4.29  at  z  =  Zj,  ?.  and 
substituting  for  Green's  functions  G,  and  G  „  frog.  Eq.  4.11  gives 

1  4 


.1 


cosh  7 '  s 
n  2 

: - r— -  pr  ( K  K  -  ?  , 

y '  s: nh  7  0  1  1 

n  n 


7  1  7.  +  P  cosh  7  '  z 

n  1  1  'n  2 


?„  s  i nh  y  '  z. 

2  n  1 


ccsh  y's2 

^7— 7777-g  [-  (KXK  +  :■  )  cosh  Yh 
n  ti 


+  P  cosh  7 ' 

1  n  1 


+  P  sinh  Y'c 
2  n  2 


t  a 


(p  -  dk2/y , 


b  -  7, 


Multiplication  of  Eqs.  1.31  by  cos  (nrn/'Ox  and  into, -ration  along  the 
interface  from  x  =  x  to  x  =  (a  -  x  }  yields 


„rr.  r  r.  ,n  1;  r. 

t  +  i  r,  f  +  ,' 

1  h  inn  '  1 1  1  1.2 


CV  +  i 


’  21  1 


(m  1  ,  X:  *  0; 


,  r.  =  n 


-  ii  gx .  a  ;n  n  ■* 


'*  '"1  i  n 


:i  =  t.-veri 


'■n/:Vjr  r.::  . 


Iqs.  U.3?  take  the  following  form  in  matrix  notation: 


To  obtain  a  iiOr.tr  i  v !  a  :  ; 

determinant  of  the  ••  - 
In  the  numerical  amlyr 
to  order  N.  If  the  f  r 
then  the  dispersion  roia 
given  by 


•  -  It  :  ‘  Cn,s,  the  infinite 

:  rix  Ir.  Zqi:. .  t.33  must  be  zero. 

infinite  determinant  in  truncated 

it  •  •••“:•.■  ::•! :it  is  defined  as  D.T(f,K), 

u 

i-  solved  numerically  for  roots  is 


=  0  .  (U.3U) 

The  discussion  on  the  r r:rer  choice  of  !i  for  reasonably  accurate 

results  and  the  numerical  data  are  presented  in  the  next  chapter. 

h  .  I4 . 2  Normal  Mag)  etj  sailor:  i  :  emulation.  The  integral  equation 

for  this  case  is  developed  by  using  the  j  roper  p.r  and  given  by 

Kqs.  4.23  ana  b.2h.  Thus ,  from  I'.q.  ‘-.2  it  follows  that 
a-x  Z  r-  -j 

!  0  {  (■  ,  I 

z;  =  ,  ,•*— -  g.’  ,z’ '  jG(x,x’  ,z,r.’  )  dz’  dx’ 


a-x  z 

— 

f  °  f  ? 

I 

f?  -  1' 

1  u 

:  v »  -t  1 1 

”;:z  I 

x  ;  z 

V.  r 

1 

0  1 

_ 

— J 

v'M  -  :)^x(xo,z’  +  K^'hr'y,.  ,z,z’)  dz' 


(u  -  1 ) <j>  /.'=  -  x  ,r  ’  •  +  _  x  ,z''it-/x,a  - 


-re  again,  the  to  it.-. 


:e-d  irr.ensiona.:  in 


r  r > x  i  .  i. o n  3  i  raj  1  i  f  i  c 

.■1  i  -  tractable. 

,  1  t::v ■  i  ig.  i'-r  r a ’ 


C'l  L  tr. 


■<  evaluate 


and  zq  =  z  =  z?  (Fig.  4.3).  The  variation  of  <j>(x,z)  in  each  layer  in 
the  z-direction  is  assumed  to  be  linear.  In  this  manner  three  functions 
each  having  one  variable,  are  used  to  approximate  <f>(x,z )  in  the  slab, 
as  follows: 

fx(x)  =  ^(x,zi)  , 

f  (x)  -  4>(x,z  ) 

o  o 

and 

f2(x)  =  ^(x,z2)  . 

The  linear  approximation  gives 


$(x,z) 


f  (x)(z  -  z)  +  f  (x)(z  -  z  ) 


z  -  z 
o  i 


5  z  =  z  =  z  (k . 36a) 

1  o 


and 


f0(x)(Z2  -  Z)  +  f2(x)(Z  -  Z0}  <  < 

d>  X  ,  z  )  2 - >  z  =  z  =  z , 

T  Z  -  Z  O  2 

2  O 


( h . 36b ) 


vliich  gives 


and 


f  ( x }  -  f  [x 

■  ,  ,  o  1  ^  < 

t>  (x,z  =  — — - - - —  ,z  =  z  =  z 

z  ’  z  -  z  1  o 


(x,z)  = 


o  1 
f  (x)  -  f  (x) 


z  -  z 
2  0 


<  < 

z  =  z  =  z 
o  0 


■;zzU’z' 


fix)  -  Of  (x)  +  f  (x) 
2  Ol 


d' 


(U.3Ta) 


( U . 3Tb ) 


(4.38) 


where  d  = 


1 


V|T; 


=  i/P  Mt:d  t  j  ;■  the  slab  thickness.  In 
.  !( .  v; ) .  the  and  third  terms  are 


the  integral 
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In  the  thin-slab  assumpt Ion ,  instead  of  a  continuous 

distribution  of  surface  charge  in  z ,  the  charge  distribution  on  both 

sides  of  the  slab  (x  and  a  -  x  )  in  each  region  z  =  z  =  z  and 

o  o  i  o 

<  < 

z  =  z  =  z  is  assumed  to  be  uniform.  This  uniform  charge  distribu- 
02 

tion  assumption  connotes  that  the  surface  charge  in  each  region 
is  equal  to  the  mean  of  its  values  at  the  edges  of  that  region. 

From  Eqs.  4.24  the  surface  charge  can  be  rewritten  as 


P 


+  (y  -  l)cj>x(x.,z)  +  K<j>(x.  ,z )  , 


( 4 . 39a ) 


where  x.  =  x  or  a  - 
i  o 

distribution  in  the 
4>(x.  ,z) 

l 


x  .  Because  of  the  uniform  surface  charge 

o 

z-directicn  in  each  layer,  it  follows  that 


2 


z  =  z 

o 


and 


<$(>:. 

l 


»z) 


f  ( X .  )  +  f  (  X  .  ) 
2  1  O  J 


Therefore,  Eq.  4.39a  for  pe(x.  ,:)  now  becomes 


T 

r;(x, ) 

+  1’^  (  X  .  )  ' 

e..'  — 

*  KxK 

[fi(xi)  + 

f0(xi)) 

+ 

|i  “  -i  ) 

2 

o 

1 

t  - 1 ) 

ff>2(Xi) 

+  f  (x.  )> 

_ 9 

1/ 

rf2(xi)  + 

ro<xih 

** 

l  2 

-  r.  r. 

1 

l  2  J 

(4.39b) 


The  main  reason  for  uniform  charge  di  stributicr.  at  the  slab  sides 

(x  =  x  and  a  -  x.  )  is  the  fact  that  <f,(x,z)  is  •  vs.  i  at 

three  values  of  z,  i  .  e .  ,  z  =  z  ,  z.r  n  n  :  .  .  he*  function  $  between 

these  values  is  unknown  and  s'  a.  i  i  t  n  •  •  ecus?  j  •  r.hou  ‘  :  involve 
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calculations  of  <Ji(x,z)  strictly  at  these  three  values  of  z.  Thus 
continuous  charge  distribution  would  not  be  a  plausible  assumption 
under  these  conditions.  Substituting  from  Eqs.  U.36  through  Eqs. 
0  39  into  the  integro-differential  equation  given  by  Eq.  U.35 
yields 


$(x,z) 


CL-“  A  Z 

f  0  f 2  fp  - 1) 

,f2(x')  -  2fo(x')  +  fi(x’)'| 

X  Jz  i 

l  d 

o  1 


G(x,x' 


)  dz '  dx '  - 


z  r 

o 


1 


(vi  -  1) 


w 


i"(x  ), 

O  O 


+  K  K 


(f, (x  )  +  f  (x  ) 
10  00 


2-  „  f— 


G(x ,xq ,z  *  )  dz!  - 


O 


(y  -  1) 


rf'ix  )  +  f'(x  ) 

2  O  O  O 


+  K  K 

1 


rf  (x  )  +  f  (x  )•* 
20  00 


G(x  ,x  ,z,z'  )  dz' 


z 

f  0 


(m  -  1) 


f  ’  (a  -  x  )  +  f '  (a  -  x  )•, 


o  o 


+  K  K 

1 


f(a-x)+f(a-x) 


o  o 


G(x,a  -  xq,  z,  z' )  dz'  + 


(y  -  1) 


f ' (a  -  x  )  +  f ' (a  -  x  ) 
2  o  oK  o 


+  K  K 

1 


f  (a  -  x  )  +  f  (a  -  x  )\ 

°  '  1 G ( x  ,a  -  x  ,z,z'  )  d z '  .  (OUo) 

o 


In  Eq.  U.UO,  the  z-inter.rals  can  be  carried  out  easily  and 


may  be  expressed  in  terms  of  on-*  or  more  of  the  followinr,  constants: 

z 


n 


ash  y  ? ' . '  d:: 


n 


1 


Substituting 
the  following 
tively : 

f\(x)  =  - 

n 


4 


cosn  v  (b  -  z  )  dz ’ 
n 
z 

o 


for  G  from  Eq.  U.ll  and  carrying  out  the  z-integrals  yields 

when  Eq.  H.40  is  evaluated  at  z  =  z  ,  z  and  z  ,  respec- 

1  o  2 


a-x 


j  a( 1  +  6  ) 

=o  on 


(Bn  +  F.n)  cosh'i  ' z  c o s  ~~  x 
'3  4  n  1  a 


nTT 

r  0 

\v  ~  1| 

X 

a 

A 

„  mr 

2v  cos  —  x 


f  (x' )  -  2f  (x’ )  +  f (x’ ) 


00 


=  -  y  — tz - - - r  [B  cosh  Y'(b  -  z  )  +  B,  cosh  y'z  ] 

L  a(l  +  6  )  l  'n  o  4  'no 

n=o  on 


r  O  ff  (x')  -  2f  (x'  )  +  f  (x') 

nit  y  -  1  2  o  i  nir 

—  x  ^ - 77 -  cos  —  x'  dx' 

a  y  d  a 

x 

o 


_  nw  r- 

2v  cos  —  x  _  (/,  ,  \ 

n  a  o  nir  n  .  I  y  —  1J  r  «,  /  \  ,  ■  r \  i 

— 7T - - — t —  cos  —  x  B,  cosh  y'z  - - -  [f  (x  )  +  f'(x  )] 

a  1  +  S  J  a  4  n  o  2  2  o  o  o 

on 


tf2(xo)  +  fo(xo)]l+  Bi  COSh  Yn(t  "  Zo)(“"i~  [fiUo] 


K  K  \  00  2v  cos  —  (a  -  x  ) 

,)]  +“T  [fi(xo)  +  fo(xo)]  -  l  n  a(  1  f-6  ~) - ~ 

n=o  on 


—  x  b"  cosh  y'z  . ■-  [fl(a  -  x  )  +  f'(a  -  x  )] 

a  4  n  o  I  2  2  o  o  o 


[  f  (a  -  x  )  +  f  (a  -  x  )  ]  +  Bl1  cosh  Y'(h  -  z  ) 

200  o  i  n  o 


v  r\  rv 

— —  ff'(a  -  x  )  +  f'(a  -  x  )1  +  — —  [  f  (a  -  x  )  +  f  (a  -  x  )' 
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f2(x) 


«  -  l 


2v 


cos  —  x  cosh  v'(b  -  z  )(RU  +  B*1) 
a  'n  21  2 


n  n , 


a(  1  +  5  ) 

n=o  on 


a-x 


o  f2(x')  -  2fQ(x ' )  +  f,(x') 


2v 


3 - dx’  -  l  IfTT 


nir 

cos  —  x 


L  a( 1  +6  )  a  Ao 

n=o  cn 


cosh  v'(b  -  z  )  cos  —  x 
rn  2  a 


(u_-  1) 


K  K 


[  f'(x  )  +  f'(x  )]  +  ~  [f  (x  ) 

2  20  OO  220 


+  f  (*  n 

o  o 


+  B 


f(u  -  1) 


K  K 


2  [f;(xo}  +  f-(xo)i  +  -f-  [fI(xQ)  +  f0(xQ)] 


2v 


)  ~ r; — v  cos  —  (a  -  x  )  cosh  v'(h  -  z  )  cos 

ia(l+6)  a  o  'n  2 

n=o  on 


n7i 


t  f 2  (a 


x  )  +  f ' (a  -  x  )] 
o  o  o 


K  K 

+  [  f 2 ( a  -  xQ)  +  fQ(a  -  xQ)] 


+  B 


(u  -  1) 

2 


[r  (a  - 


Xo} 


+  f' 


Xo)] 


XiK 


[fx(a  - 


xo} 


fo(a  - 


Xo^ 


(i* .  hlo ) 


Equations  describe  a  set  of  three  coupled  equations  in  terms  of 

f  ,  f  and  f  .  By  introducing-  the  following  functions,  this  set  can 
10  2 

be  reduced  to  a  set  of  two  coupled  equations  which  is  more  attractive: 

G  (x)  =  f  (x)  -  2f  (x)  +  f  (x)  ( b .  l+2a ) 

V  2  0  1 

and 

G  (x)  =  f  (x)  +  2f  (x)  +  f  (x)  .  (h.k2b) 

S  2  0  1 

Forming  these  functions  gives 


r2(x)  -  2fo(x)  +  fx(x)  = 


1  ~~rr~zh — y cos  ir  xf(B,  +  )  cosh  Y' 

t  a(  1  +  6  )  a  12  r 


(b  -  z  )  +  (B"  +  B")  cosh  -  2b"  cosh  y^b  -  zq)  -  2B^  cosh  y^ 


[  f  (x'  )  -  2f  (x'  )  +  f  (x'  )  ]  COS  —  x'  dx'  +  cos  ~  x 

U  o  2  O  1  3-  dL. 

Jx  d  J 

_  o 

K  K 

(-  [f2<*o' ♦  2fouo»  *  w]  -  -r  (f;uo> +  2f;(xo» 


,  .  ,  I  nn (a  -  x  )  ,  ,  \ 

+  fj(x  )]  +  COS  - o_  ,(W  -  u  [f  (a  _  x  )  +  2f  (a  -  x  ) 

1  °  ;  a  i*  L  2 1  o  o  o 


fi(a-xo)]  +  ~jr  -  x0}  +  2fi(a  -  xo}  +  fi(a  "  xo)]J  {U'k3&) 


f  u)  +  2f  (x)  +  fjU)  =  -  i  TfrrV^y cos  r  xC(Bi  +  b2} 

n=o  on 

•  cosh  Y'(b  -  z,)  +  (B*  +  Ba)  cosh  y'z  +  2Ba  cosh  y'(b  -  z  ) 
z  3h  n  a  A  'w/ 


a-x 

n  I  °  (u  -  I)  r  „ 

2B  cosh  y'z]  ~  [f, 

4  n  °  Jx  V 

o 


[f2(x>)  -  2fQ(x’)  +  fj(x’)l 


an  ,  ,  ,  .  „  1111  v  in  -  ^ 

cos  —  x'  dx'  +  cos  — -  x  -  r 
a  a  o  u 


[f„(x  )  +  2f  (x  )  +  f  ( x  )  ] 
1  2  '  o  O  O  10 


nlT(a  -  X0}f(u  -  1) 


_  -4-  [  f '  (x  )  +  2f '  (x  )  +  f '  (x  )  1  +  cos 
h  1  2  o  o  o  1  o 


[f2(a  "  X0}  +  2fOU  '  X0}  +  fl(a  “  XG)]  +  "V  [f2(a  ‘  X0} 


+  2f'(a  -  x  )  +  f'(a  -  x  )]  .  (U.k3b) 


By  defining,  the  following  constants,  Eqs.  U .  U  3  simplifies  greatly: 


v  [ (fin  +  Bn)  cosh  v'(b  - 
n  1  2  n 


)  +  (  B1'  +  l?n )  cosh  y  '  z, 
3  4  n  i 


+  2B  cosh  y ' ( b  - 
l  n 


)  +  2B‘  cosh  y'z 
4  no 


W 


v  [  (B.n  +  B’1 )  cosh  v '  (b  -  z  )  +  ( Bn  +  i>n )  cosh  y'z 
n  l  2  n  2  3  4  ni 


-  2BR  cosh  Y'(b  -  z  )  -  2fcn  cosh  y'z  ] 
l  n  o  4  'no 


and 


a-x 
r  o 


V  V 


0  (x)  cos  —  x  dx 
a 


nn 

cos  -  X 

a  o 


h  -  1 


K  K  dG  (x  ) 


,  G  (x  )  +  — r—  — 

4  S  o  4  QX 


nn  , 

4  cos  - —  (a-x  ) 
a  o 


,  K  h  dG  ( a  -  x  ) 

y  -  1  ,,  ,•  \  .  i  s  o 

: -  G  (a  -  x  )  +  —j—  - - - 

4  s  o4  dx 


Thus  Eqs.  4.43  can  be  concisely  stated  as 


G  (x) 
v 


-  I  r 


a(  1  +  6  ) 

n=o  on 


nn 

COS  -  X- 

a 


(4.44a) 


and 


.(x)  =  -  l 


roTTT  cos  r 

n=o  on 


(4.44b) 


The  term  Qn  in  Eqs.  4.44  involves  the  first-order  derivative  of 
G  (x).  The  function  G  (x)  [obtained  from  f  (x),  f  (x)  and  f  (x)] 

s  S  1  O  2 

is  defined  to  be  nonzero  in  the  ferrite  slab  (x  =  x  =  a  -  x  )  and 

o  o 


zero  everywhere  outside*  which  nieuns  ‘that,  :.t  is  discontinuous  at  x  — 
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and  x  =  a  -  and  thus  its  derivatives  in  the  x-direction  at  the  slab 

edges  (x^ora  -  x  )  are  undefined.  This  problem  creates  difficulty  in 

the  evaluation  of  Qn .  To  overcome  this  problem  of  discontinuity,  the 

function  Gg(x)  is  defined  only  in  the  range  xq  =  x  =  a  -  xq.  In  this 

manner  the  function  G  (x)  becomes  differentiable  at  x  =  x  or  a  -  x 

s  o  o 

and  its  approximate  finite  series  expansion  can  be  written  as  follows: 


G  (x, 
s 


=  p 


N 

o 

l 

a  = i 


cos 


l  n 


2x 


(x  -  x  ) 
o 


N 

o 

-  I 

1=1 


sm 


2x 


(x  -  x  ) 
o 


where  p^'s 

dG  ( x) 
s 

dx 


x  =  x  =  a  -  x  ,  (U.U5a) 

o  o 

and  q^'s  are  arbitrary  constants.  Upon  differentiation. 


I 


*=1 


S,  u  .  i  n 

- - —  sin  - - — 

a  -  2x  a  -  2x 

o  o 


(x-xo) 


N 

,  r  U 

L  %  a  -  2x 
1=1  o 


cos  — — — —  (x-x)  x  =  x  =  a  -  x  ,  (U.l*5b) 

a  -  2x  o  o  o 


where  H  is  a  very  large  integer  number.  Therefore,  when  the  series 

expansions  for  G  (x)  and  [  dG  (x)/dx]  are  used,  the  expression  for  0° 
s  s 

becomes 


a-x 


n  u  -  1 

Q  =  ,2 

pd 


,  /  i  rrn  .  nn 

'  i  [  X  cos  -  X  iX  -  COS  —  X 

v  h  a  o 


P  -  1 


Po  +  l 
{  1=1 


KK  No 

l  Y  iv 

.1 

Tin  /  \ 

+  cos  —  ( a  -  x  ) 

*  i  0 

p  - 1 

*  I  (-  d'p, 

+  a  l  a  -  MX 

*  »' 

y.=l 

__ 

l  «=1 

'•  ,  ,  >1ti 

L  ^  - - -  G;, 


a  -  c1  x  i. 


.  (U.U6) 


To  obtain  Use  dlspers 
Mu  It  i  t  1  i  cat  i  t  p 


ioi  rt-lati  n,  the  following  procedure  is 
f  iv  cor  (itn/a)x  and 


undertaken : 


sin  (mir/a)x  and  integration  in  x  from  x  =  x  to  a  -  x  yields  the 

o  o 

following  linearly  independent  equations: 


where 


C  +  £  a 

V  u  TT!  >1 

n=o 


(U.UTa) 


-.It.  r  n,  n 

n  o  +  l  a  4  o  =  0 


(U.l+Jb) 


sm  +  l  e  q¥  =  o  , 

C  L  rr.n 
n=o 


( h . hTc ) 


mt  „  /  %  . 

cos  —  x. i  x)  ox  , 
a  v 


III:!  /  \  , 

ccs  —  x.'j  J.x  j  dx  , 
a  o 


•  mu  ,,  fv)  dx 
sm  —  x-d  y 


d  rmt  rni 

— tt - - t  cos  —  x  cos  —  x  dx 

all  +  6  )  a  a 

on  x 

o 


d  .  rniT  n-rr  , 

—ro - — : - 1  sm  —  x  cos  —  x  dx  , 

all  +  6  )  a  a 

on  ^  x 

o 


n  =  0  1  2  N 


ITi 


I’he  term  Q  from  Eq.  U.Uf  is  expressed  in  terms  of  constant  coefi'i- 


.  but  as  s‘f-n  from  th»-  !  : .me  :  inoarlv  inlet endent 


cients  t>  and  o 
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equations  (Eqs .  4.1+7),  it  is  best  to  describe  Qn  in  terms  of  c‘",  C™ 

v  s 

and  Sm. 
s 

To  find  p  ' and  q  '  s  in  terms  C™  and  S^,  Eq.  4.45a  for  f!_(x) 

X/  X/  S  S  o 

is  multiplied  by  cos  (rrm/a)x  and  sin  (rrnr/a)x,  respectively.  The 

resulting  equations  are  then  integrated  from  x=x  to  x  =  a  -  x  . 

o  o 

Thus  the  following  two  equations  are  obtained: 


(cc)mp  +  (cc)m  +  (CC)mp  + 

0  0  1  2  2 


.  +  (SC)mq  +  (SC)mq  +  . 
11  2  2 


( 4 . 48a ) 


(CS)^Pc  +  (CS)“Pl  +  (CS)”p2+...  +  (SS)“qi  +  (SS)^q2  + 


(4.48b) 


where 


nit  ,  ,  mu 

)  =  cos  — - - —  (x  -  x  )  cos  —  x  dx  , 

n  a  -  2x  o  a 

l  x  c 

o 


)  =  cos  — — - —  x  -  x  sin  —  x  dx  , 

n  a  -  2x  o  a 

J  X  o 

o 


vm  riTi  /  v  imr 

)  =  sin  - - —  (x  -  x  )  cos  —  x  dx  , 

n  a  -  2x  o  a 

;  x  o 

o 


xm  mr  (  \  in.! 

)  =  sm  - r —  lx  -  x)  r-in  —  x  dx  , 

n  a  -  2x  o  a 

J  x  o 

o 


n  =  0,1,2,.  .  .  ,  N 


m  =  0,l,2,...,i; 


Writing 


in  matrix  format 


v  -i  e  l  o  c 


C'' 


s 


0 

*  o 


s 


r> 


1 


H  ]  • 


coefficient  matrix  given  b; 
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Premultiplying  Eq.  4.49  by  [H]-1  and  observing  [ H ] [ H ] — 1  =  1 

yields 


(4.50) 


Substitution  of  p  's  and  a,  's  in  terms  of  Cm,s  and  Sm,s  in  Qn 

of  Eqs.  4.47  will  produce  a  system  of  linear  equations  in  C™,  and 

Sm.  To  obtain  a  nontrivial  solution  for  this  system  of  linear  equa- 
s 

tions  requires  that  the  large  determinant  (Nq  x  N  )  of  the  coeffi¬ 
cient  matrix  be  set  to  2ero. 

In  practice  the  actual  size  of  the  matrix  needed  to  obtain 

a  reasonably  accurate  answer  is  smaller  than  (N  x  N  )  and  it  depends 

00 

on  the  particular  mode  under  consideration.  In  Chapter  V  Eqs.  4.47 
are  discussed  in  detail  for  the  first-order  mode  (m  =  l)  and  a 
computer  program  based  on  a  truncated  matrix  is  written  to  provide 


numerical  insight  into  the  problem. 


CHAPTER  V.  COMPUTER  SIMULATION  AND  RESULTS 


3.1  Introduction 

The  mathematical  formulations  presented  in  Chapters  III  and 
IV  were  programmed  on  a  computer  and  the  effect  of  different  param¬ 
eters,  such  as  slab  position,  width  and  thickness,  on  the  dispersion 
characteristics  and  group  time  delay  was  investigated.  Most  of  the 
numerical  results  were  obtained  for  a  fixed  dc  field  value  of 
H^c  =  1800  Oe  and  a  YIG  saturation  magnetization  of  M^  =  1750  Oe. 

The  tunable  properties  for  various  magnetic  dc  fields  were  also 
investigated  and  the  results  were  presented. 

As  seen  earlier,  calculation  of  the  dispersion  relation  for 
the  purpose  of  root  finding  involves  evaluation  of  large  matrices. 

In  Section  5*2  algorithms  to  find  the  determinant  and  inverse  of  a 
matrix  are  presented.  This  procedure  was  programmed  for  use  in  all 
of  the  determinant  calculationo. 

Root  finding  of  the  dispersion  relation  is  done  by  the 
Newton-Raphson  method  which  is  described  in  Section  5-3.  A  flow 
chart  is  also  provided  for  use  in  programming.  With  the  aid  of 
this  algorithm,  determinant  roots  of  the  dispersion  relation  are 
found  through  several  iterations. 

In  Section  5*^,  the  dispersion  plots  corresponding  to  the 
mode  analysis  for  both  parallel  and  normal  magnetization  are 
presented  and  several  of  its  ramifications  and  consequences  are 


described. 


Results  for  the  integral  equation  method  are  presented  in 
Section  5*5-  The  dispersion  plots  and  group  time  delay  are  pre¬ 
sented  and  discussed  for  both  parallel  and  normal  magnetizations. 
For  the  normal  magnetization  case  the  analysis  for  the  first-  and 
zeroth-order  modes  and  its  programming  are  presented  in  detail. 


is  now  the  largest  of  the 


set  { jaK).|  ,  la^  \  ,..., 


3.  Step  3.  Test  for  near-zero  maximum  element  (singular  or 


near-singular  matrix).  If  |a^|  =  e  go  to  error  exit.  If  la^l  >  e 
continue  to  Step  4. 

4.  Step  4.  Perform  stage  K  of  reduction  process: 


aKK’ 


a.  Row  K  =  (Row  K/a^^)  and  D  =  aT  «D.  If  div  = 

Kiv  KK 

then  aRj  =  a^/div  and  b^  =  b  /div  (j  =  1,  n). 

b.  Row  i  =  Row  i  -  a...  Row  K.  If  mult  =  a.„,  then 

lK  xK 

a .  .  =  a .  .  -  mult  •  a„  .  and  b .  .  =  b .  .  -  mult  »bt, , . 

Ij  ij  aj  ij  ij  KJ 

5.  Step  5.  Test  stage  counter  K.  If  K  <  n,  set  K  =  K  +  1 
and  return  to  Step  2.  If  K  =  n  continue  to  Step  6. 

6.  Step  6.  Write  output  B  =  A-1  and  D  =  [ A I . 


3.3  Newton-Raphson  Method 

The  equation  to  be  solved  numerically  is  the  dispersion  rela¬ 
tion,  a  function  of  frequency  (f)  and  wave  number  (K) ,  which  can  be 
written  as 

D(f,K)  =  0  ,  (5.1) 


where  D(f,K)  represents  the  determinant  of  the  coefficient  matrix 
involved  in  the  system  of  linear  equations.  Equation  5.1>  in 
general,  is  a  nonlinear  function  of  f  and  K  and  can  be  quite  compli¬ 
cated  if  the  size  of  the  matrix  is  large. 

With  the  aid  of  the  Newton-Raphson  method,  Eq.  5.1  is 
solved  numerically  for  roots  K  (at  a  known  frequency  f^).  The  follow¬ 
ing  algorithm  details  the  exact  steps  used  in  programming  Eq.  5>1 
in  order  to  find  its  roots: 

1.  Step  1.  Input  and  Definition.  Read  f  =  the  frequency  of 
operation,  K  =  the  initial  approximation  of  the  root  of  D(fi,K)  =  0, 


e  =  the  convergence  term  and  N  =  the  maximum  number  of  iterations. 

2.  Step  2.  Initialization.  Set  iteration  counter  i  =  1.  Set 

the  correction  term  A  =  C  (C  is  an  arbitrary  large  positive  number), 

111 

3.  Step  3.  Compute  successive  approximation  of  root  using  the 
Nevton-Raphson  iterative  formula: 


i+i 


K. 

l 


D(fx,K.) 
D'Cf^K.)  * 


Compute  the  magnitude  of  the  correction  term  in  the  current  iteration 


Ai+: 


K. 


i+i 


-Kil 


4.  Step  4.  Test  for  convergence  or  failure  to  converge. 

a.  If  Ai+i  =  e  and  lD(f1»Ki+1)|  =  e *  g°  to  Step  5.  If  not 

continue. 

b.  If  >  A^,  select  new  and  return  to  Step  2.  if 

A .  =  A .  continue, 

i+i  i 

c.  If  i  *  N,  set  i  =  i  +  1  and  return  to  Step  3.  If  i  >  N, 
select  new  K  and  return  to  Step  2. 

5 .  Step  5.  Output  root  Kq;  set  Kq  =  K^+i .  Write  Kq.  Steps 
1  to  5  are  summarized  in  the  flow  chart  shown  in  Fig.  5.1. 


5.4  Mode  Analysis  Computer  Simulation 

In  Chapter  III  magnetostatic-wave  propagation  in  a  waveguide 
partially  filled  with  a  YIG  slab  was  analyzed.  The  YIG  slab  was  in 
contact  with  the  sidewalls  of  the  waveguide.  The  dispersion  relations 
for  the  two  principal  directions  of  magnetization  were  formulated. 

For  parallel  magnetization ,  an  exact  equation  for  the  disper¬ 


sion  relation  was  derived  for  a  YIG  slab  of  arbitrary  thickness. 
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For  normal  magnetization,  there  is  no  exact  solution  possible. 

However,  the  dispersion  relation  in  this  case  is  obtained  in  terms 
of  an  infinite  determinant  that  should  be  truncated  properly.  The 
order  of  truncation  depends  primarily  on  the  particular  mode  (m) 
under  consideration.  Numerical  root  findings  of  either  the  exact 
equation  or  the  truncated  determinant  on  a  computer  are  facilitated 
with  the  use  of  the  Newton-Raphson  iterative  algorithm  given  in 
Section  5.3.  The  corresponding  numerical  data  for  the  w-k  diagram 
and  group  time  delay  vs.  frequency  for  parallel  and  normal  magneti¬ 
zation  directions  are  presented  and  discussed  in  the  following 
sections . 

5.U.1  Parallel  Magnetization.  From  Eq.  3. IT,  several  effects 
in  terms  of  different  device  dimensions  may  be  studied.  The  wave¬ 
guide  dimensions  are  assumed  to  be  2.0  x  1.0  cm  and  the  slab  thick¬ 
ness  is  0.1  cm.  The  dispersion  characteristics  are  plotted  in  Fig. 
5.2.  In  this  figure,  as  the  YIG  slab  is  lowered  from  the  top  to  the 
bottom  of  the  waveguide,  the  dispersion  characteristics  are  restricted 
tc  lower  propagation  bandwidths.  The  middle  positions  of  the  slab 
in  the  upper  half  of  the  guide  in  some  frequency  ranges  exhibit  nega¬ 
tive  slope  corresponding  to  negative  group  velocity  which  means 
energy  propagates  in  the  opposite  direction  to  wave  propagation.  The 
corresponding  group  time  delay  in  ns/cm  defined  by  the  relation 
=  (3u)/3K)-1  is  shown  in  Fig.  5.3.  From  this  figure,  it  can  be 
seen  that  as  the  slab  separation  from  the  bottom  surface  is  reduced, 
the  time  delay  increases  while  the  bandwidth  decreases.  It.  should 
be  noted  that  these  figures  are  drawn  for  the  first  mode  and  only 
for  the  forward  propagating  wave. 


(a)  zq  =  0.95  cm,  (b)  zq  =  0.8  cm, 

(c)  z  =0.7  cm,  and  (d)  z  =  0.05  cm 
o  o 


FIG.  5-2  EFFECT  OF  SLAB  POSITION  ON  THE  DISPERSION 


CHARACTERISTICS. 


FREQUENCY,  GHz 


EFFECT  OF  SLAB  POSITION  ON  GROUP  TIME  DELAY  PER  UNIT 


Figures  5.^  and  5-5  show  the  effect  of  slab  thickness  on  the 


dispersion  characteristics  and  group  time  delay  per  unit  length. 

As  the  slab  thickness  is  reduced  the  wave  propagate  in  a  smaller 
bandwidth  but  with  a  larger  group  delay  (first  mode). 

Figure  5.6  shows  the  first-  and  second-order  modes  for 
different  slab  positions.  As  can  be  seen,  higher-order  modes  have 
cut-off  points  on  the  K-axis  and  as  the  mode  number  increases,  the 
cut-off  point  also  increases.  Modes  higher  than  the  second  are  not 
town  in  this  figure. 

Tunable  properties  are  also  investigated  by  varying  the  magnetic 
bias  field.  Figure  5-7  shows  the  effect  of  magnetic  bias  field  on  the 
dispersion  curves.  As  can  be  seen,  the  dispersion  curves  move  up  or 
down  the  w-k  plane  by  varying  H^.  Figure  5.8  shows  the  corresponding 
effect  of  bias  field  on  the  time-delay  characteristics. 

5.^.2  Normal  Magnetization.  The  zeroth  order  mode  is  already 
discussed  in  Section  3.5.  To  obtain  a  nontrivial  solution  for  higher 
order  modes  from  the  system  of  linear  equations  given  by  Eq.  3.22, 
the  determinant  of  the  coefficient  matrix  which  is  infinite  in  size 
must  be  set  to  zero.  However,  for  practical  purposes,  the  matrix  was 
properly  truncated  for  best  accuracy.  The  truncation  cut-off  point  of 
the  matrix  depends  on  the  mode  of  propagation.  For  example,  for  the 
first-  and  second-order  modes,  matrices  of  orders  up  to  12  x  12  were 
studied  and  minimum  matrix  size  was  found  to  be  U  x  U  and  6x6,  respec¬ 
tively.  For  higher-order  modes,  larger  matrices  must  be  considered. 

A  computer  program  was  written,  using  the  determinant  algorithm 
presented  in  Section  5*2,  to  find  the  determinant  of  the  truncated 
coefficient  matrix.  Roots  of  this  equation  are  found  by  using  the 
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Newton-Raphson  method  as  described  in  Section  5.3.  Through 
numerical  analysis  it  is  found  that  the  magnetostatic-wave 
propagation  is  symmetrical  in  the  guide  cross  section  with  respect 
to  the  slab  position.  This  observation  was  made  possible  by 
noticing  that  the  dispersion  curves  for  slabs  positioned  in  the  upper 
half  of  the  guide  are  the  same  as  for  slabs  positioned 
symmetrically  in  the  lower  half  of  the  guide.  Unlike  the  parallel 
magnetization  case,  the  wave  propagation  is  reciprocal  for  K  and  -K. 

Figure  5-9  shows  the  effect  of  lowering  the  slab  position. 

The  propagation  bandwidth  is  reduced  as  can  be  seen  from  this  figure. 
Besides  the  zeroth-order  mode,  two  higher-order  modes  are  also 
shown.  These  higher-order  modes  exist  due  to  the  finite  width  of 
the  slab.  It  is  noted  that  there  is  a  frequency  cut  off  for 
higher-order  modes  which  increases  as  the  mode  number  increases. 

In  Fig.  5*10,  the  group  time  delays  for  different  modes  are 
plotted  when  the  slab  is  placed  against  the  top  or  bottom  of  the 
guide.  This  figure  shows  the  increase  in  time  delay  as  the  mode 
number  increases.  Higher-order  modes  have  higher  frequency  cut¬ 
off  points. 

Figure  5.11  shows  the  effect  of  the  slab  position  on  the 
group  time  delay  per  unit  length.  The  slab  is  placed  in  two 
positions,  at  the  center  of  the  guide  and  at  the  bottom  or  top  of  the 
guide.  The  time  delay  increases  as  the  slab  is  moved  toward  the 
center  of  the  guide. 

Figure  5-12  shows  the  effect  of  normal  magnetic  field  for 
several  bias  field  values.  This  effect  on  the  time-delay  charac¬ 
teristics  is  shown  in  Fig.  5.13. 
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?IG.  5-13  GROUP  TIME  DELAY  CHARACTERISTICS. 


.5  Integral  Equation  Computer  Simul at i on 


In  Chapter  IV  the  integral  equation  method  using  Green's 
function  was  applied  to  the  problem  of  magnetostatic-wave  propaga¬ 
tion  in  a  YIG  slab  of  finite  width  in  a  waveguide.  The  dispersion 
relations  for  two  principal  directions  of  magnetization  were  derived. 
These  dispersion  relations  are  in  terms  of  infinitely  large  deter¬ 
minants  as  given  by  Eqs.  i*.33  and  h.kj. 

To  be  able  to  solve  these  determinant  equations  on  the  computer, 
they  must  be  properly  truncated  at  some  cut-off  point.  This  trun¬ 
cation  problem  was  studied  on  the  computer  by  the  method  of  trial 
and  error.  As  a  result,  the  trend  of  choosing  the  proper  matrix 
size  was  obtained.  Thus,  it  is  found  that  the  order  of  truncation 
depends  on  the  particular  mode  under  consideration  (m)  and  the  size 
of  the  air  gap  existing  on  either  side  of  the  YIG  slab  (xo). 

Several  theoretical  results  showing  this  correlation  are  presented. 

Since  parallel  and  normal  magnetization  each  involve  separate 
formulations,  each  case  is  discussed  and  programmed  separately  and 
the  results  are  presented. 

3.3.1  Parallel  Magnetization  Results.  To  obtain  a  nontrivial 
unique  solution  for  Cn,s  in  Eq.  U.33,  the  infinite  determinant  of 
the  coefficient  matrix  is  set  to  zero  in  theory.  However  in  practice, 
the  size  of  the  matrix  is  reduced  by  finding  a  proper  cut-off  point 
(N).  When  the  determinant  of  the  cut  off  (N  x  N)  matrix  is  Djj(f,K), 
the  equation  to  be  solved  can  be  written  as  D  (f,K)  =  0. 

With  the  aid  of  a  computer  program,  results  were  obtained  for 
a  YIG  slab  against  the  upper  surface  of  the  waveguide  with  no  wall 
gap  (i.e.,  xq  =  0.0).  The  results  are  presented  in  Table  5.1  for 


Table  5.1 

Comparison  of  Results 
N  =  8 


Relative  Error 


Frequency 
f  (GHz) 

Exact  Root 

K  ( cm~ 1 ) 

0 

Approximate  Root 

K'  (cm-1) 

0 

[ (K*  -  K  )/K  ] 
0  00 

( Percent ) 

7.5 

-  1.145 

-  1.147 

0.14 

7-9 

-  2.050 

-  2.053 

0.15 

8.9 

-  5-379 

-  5.382 

0.05 

9.5 

-  9.7^8 

-  9.760 

0.12 

9.9 

-  22.843 

-  23.011 

0.73 

a  =  2  cm 
b  =  1.0  cm 
t  =  O.t  cm 
Hdc=  1800  Oe 


N  =  8  and  are  compared  with  the  exact  results  from  the  earlier  mode 
analysis.  For  all  practical  purposes,  it  can  be  seen  from  Table  5.1 
that  the  percent  error  in  all  cases  is  less  than  one  percent  and 
therefore  the  integral  equation  does  offer  a  very  reliable  and 
accurate  method  to  calculate  dispersion  curves,  as  far  as  thin  films 
are  concerned. 

The  effect  of  increasing  the  cut-off  point  (N)  (from  2  to  ^0) 
on  the  dispersion  curves  is  shown  in  Fig.  5.lU.  From  this  figure  it 
can  be  seen  that  for  smaller  wall  gaps,  a  smaller  matrix  size  is 
needed  to  produce  an  accurate  result. 

Figure  5.15  shows  the  relationship  of  the  normalized  wall  gap 
(2xQ/a)  and  the  cut-off  point  (N).  Roughly,  there  is  an  exponential 
increase  in  matrix  size  as  the  normalized  wall  gap  (2xQ/a)  increases. 
This  indicates  that  for  narrower  slabs  more  terms  are  needed  to  yield 
an  accurate  dispersion  relation  which  can  be  verified  mathematically. 

The  combined  effect  of  the  position  and  width  of  the  YIG  slab 
is  shown  in  Fig.  5 • l6 .  The  time  delay  is  plotted  in  Fig.  5*17.  It 
can  be  seen  that  as  the  slab  width  decreases  the  delay  time  increases 
and  the  dispersion  curves  bandwidth  shifts  downward,  while  as  the 
slab  position  is  lowered,  the  delay  time  increases  and  the  dispersion 
curves  are  compressed  with  smaller  bandvxdths.  This  means  that, 
roughly  speaking,  the  position  of  the  slab  controls  the  bandwidth  and 
its  width  controls  the  center  frequency  of  the  device.  Figure  5-18 
shows  the  effect  of  increasing  the  normalized  air  gap  (2x  /a)  on  the 
dispersion  relation.  It  is  noted  that  the  dispersion  curves  shift 
downward  as  the  air  gap  increases.  Figure  '.19  plots  the  wave  number 
vs.  the  normalized  air  gap  and  the  frequency  is  used  as  a  parameter. 
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FIG.  5.15  RELATIONSHIP  OF  THE  WALL  GAP  (xq)  AND  THE  CUT-OFF  POINT  (N). 
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FIG.  5.19  WAVE  NUMBER  VS.  THE  NORMALIZED  AIR  GAP. 


From  this  figure  it  c*u.  he  seen  *  hat  an  increase  in  frequency  leads 
to  a  higher  wave  nutter  \i.e.,  smaller  wavelength)  at  a  particular 
x  .  This  figure  also  shows  *.  hat  at  a  certain  frequency,  wavelength 
does  not  change  much,  as  the  air  gai  (x.)  is  increased.  This  fact  can 
be  used  in  transducer  design  ?••  r  excitation  of  a  certain  wavelength 
at  a  particular  frequency. 

5.5- 2  Formal  Magnetization  Results.  In  Section  U.U.2  the 
basic  formulation  for  magnetostatic-wave  propagation  for  a  normal 
bias  field  was  derived  and  was  given  by  Eiqs .  U.U7.  In  this  section 
only  the  zeroth-  and  first-order  modes  are  considered.  Equations 
h . U7  are  used  to  derive  the  dispersion  relations  in  terms  of  the 
determinant  of  a  (5  x  c)  matrix.  The  details  of  this  derivation  are 
given  next  for  more  clarity. 

For  the  first-  and  zercth-ord er  mode  (m  =  0,  1  and  n  =  0,  l) 
Eqs.  ^.1*7  become 


C°  +  a  Q°W°  +  a  ij’v1  =  0 

V  CO  '  Cl  ' 


(5.2a) 


C1  +  a  Q°W°  +  a  Q1 W1  =  0 

v  10  11 


(5.2b) 


C°  +  a  CV  +  a  Ou1  =  0 

s  00  '  01  ' 


(5-2c) 


and 


C1  +  a  Q°U°  +  a  Q 1  (J 3  =  0 

s  10  11 


C1  +  [•'  Q°U°  +6  Oil1  =  0 

s  10  '  11' 


(5. 2d) 


( 5 • 2e  ) 


where 


^  Jr  pi 


C1  -  cos 


JL  „  IP  -  -l 


yd‘  v  a  Xo  '  4  ^  '^o  T  yi '  T  4(a  -  2x  )  qi 

v  o 


(p„  +  P,  )  + 


+  cos  —  (a  -  x  ) 
a  o 


oo 


01 


10 


and 


io 


1 1 


P.  -  1 


K2K  u 


[  4  (Po  +  Pi }  4(a  -  2x  )  qi 

'  o 


,  (5.3b) 


1  -  2x  /a  , 
o 


sin  —  (a  -  x  ) 
a  o 


sin  —  x 
a  o 


01 

2 


2x  . 

a  =  1  -  —  -  +  i- 

li  a  2tt 


X  > 

sin  2n 

M 

*x  (  x 

o  ,  o 

cos  -  -  COS  It  1  -  — 

a  I  a 


1_ 

2  it 


cos  2tr  —  -  cos  2tt 
a 


x  t 

1  - 
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As  pointed  out  in  Section  4.4.2,  the  term  Qn  as  given  by  Eq. 


4.44  is  in  terms  of  and  .  To  be  able  to  solve  Eqs.  5.2  for  the 
dispersion  relation,  Q°  and  Q1  must  be  expressed  in  terms  of  more 


familiar  constants  C°,  C1  and  S1.  From  Section  4.4.2,  for  m  =  0,  1 

S  S  8  ’  ’ 


and  n  =  0,  1  it  follows  that 
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(SS)* 


a-x 

o 

sin  - - —  (x  -  x  )  sin  —  x  dx 

a  -  2x  o  a 

x  o 
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1 

2 


a  -  2x 


TT_ 
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sm 


TTX 


a  -  2x 


a  -  2x 


1 

2 


a 


■ -  +  — 

2x  a 


sm 


7IX 


a  -  2x 


a  -  2x 


n 

+  — 

a 


a-x 


To  express  p.  and  q  in  terms  of  Cm  and  Sm,  the  inverse  of  matrix  [H], 
3C  S  S 

i.e.,  [ H ]  1  must  he  known.  Matrix  inversion  is  done  with  the  aid  of 
a  computer,  so  the  results  of  this  inversion  are  used  in  the  follow¬ 
ing  equations.  To  proceed  with  the  analysis,  assume  [ H ]  1  is  known 
and  let 
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hil 

h;3 

h21 

h' 

2  2 

h' 

2  3 

• 

(5.6) 

h' 

3  1 

h’ 

32 

h’ 

3  3  _ 

5.3  for  Q1 

and  managing  the 

terms 

,'ields 


■~r  c1  +  D  C°  +  P  C1  +  D  s1 
pd^  v  os  is  2  s 
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where 
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Thus  with  the  help  of  Eqs.  5-34  and  5.7  for  Q°  and  Q1 ,  Eqs.  5-2  in 
matrix  notation  become 


[  M(  f ,K  )  ] 


=  0  , 


(5.8) 
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Requiring  a  nontrivial  unique  solution  yields  the  dispersion  relation. 
This  dispersion  relation  is  obtained  by  setting  the  determinant  of 
[M(f,K)]  to  zero. 

For  the  special  case  of  the  zeroth-order  mode,  the  dispersion 
relation  is  given  by  (n  =  0,  n  =  0): 

1  +  FaooW°  =  0  .  (5.10) 

To  find  the  dispersion  relation  for  the  first-order  mode,  the 
following  equation  must  be  solved: 

| m( f ,K) |  =  o  .  (5.11) 

Separate  computer  programs  were  written  for  each  mode  to  calculate 
the  roots  of  Eqs.  5*10  and  5.11* 

Figure  5.20  shows  the  effect  of  slab  position  in  the  waveguide 
on  the  dispersion  characteristics.  From  this  figure  it  can  be  seen 
that  the  effect  of  slab  position  on  the  dispersion  curve  becomes 
pronounced  at  the  higher  frequencies  in  the  propagation  band. 

Although  the  characteristics  all  converge  at  the  lower  end  of  the 
propagation  band,  their  slopes  are  different.  This  leads  to 
different  group  time  delays  as  can  be  seen  in  Fig.  5.21  •  This  figure 
shows  the  group  time  delay  corresponding  to  Fig.  L.?0.  As  noted  in 
Fig.  5-21,  as  the  slab  is  placed  toward  the  center  of  the  guide,  the 
group  time  delay  increases  while  the  propagation  bandwidth  decreases. 

Width  effects  on  the  device  performance  was  also  studied  and  Lh 
results  are  shown  in  Figs.  5-22  and  5.23.  In  Fig.  5.22,  it  can  be 
seen  that  as  the  normalized  air  gap  increases  the  propagation  band¬ 
width  decreases  and  the  curves  flatten  'ut  ns  they  shift  toward 


FIG.  5.20  DISPERSION  CURVES  FO! 
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5.23  TIME  DELAY 


DIFFERENT  S 


higher  frequencies.  Figure  5.23  shows  the  corresponding  group  time 
delay  vs.  frequency.  From  this  figure  it  can  be  seen  that  as  the 
slab  width  decreases  (or  the  air  gap  increases)  the  group  time  delay 
increases  toward  higher  values  with  smaller  bandwidths  as  noted 
earlier.  The  group  time  delay  at  smaller  slab  widths  remains 
constant  in  a  larger  bandwidth  and  also  has  a  higher  value.  This 
property  can  be  used  effectively  in  device  design  to  obtain  a  con¬ 
stant,  high  group  delay  per  unit  length  in  a  desired  frequency 
band.  Figure  5.2k  plots  wave  number  K  vs.  the  normalized  air  gap 
(2xQ/a).  In  this  figure,  the  information  of  Fig.  5.22  is  rearranged 
in  a  different  fashion.  It  can  be  seen  that  the  wave  propagation 
at  small  slab  widths  (or  large  air  gap)  is  possible  only  at  higher 
frequencies  with  smaller  wavelengths  (or  higher  K).  Once  the  slab 
width  is  chosen,  Fig.  5-2U  shows  the  frequency  at  which  the  device 
must  be  operated  to  obtain  a  certain  wavelength,  and  vice  versa. 

Thickness  effects  are  shown  in  Fig.  5.25.  This  figure  shows 
that  as  the  normalized  thickness  (t/b)  increases  the  dispersion 
characteristics  shift  upward  toward  higher  frequencies  with  lower 
propagation  bandwidths.  The  corresponding  group  time  delay  per  unit 
length  is  shown  in  Fig.  5.26.  In  this  figure  it  is  noted  that  as  the 
normalized  slab  thickness  increases  from  zero,  the  group  time  delay 
decreases  and  achieves  its  minimum  value  at  t/b  =  0.5  (i.e.,  when 
50  percent  of  the  guide  is  filled  with  YIG).  From  this  point  on,  as 
the  slab  thickness  increases,  the  time  delay  increases  rapidly  to 
high  values.  Although  the  time  delay  value  is  very  high  for 
0.5  <  t/b  <  1,  the  available  bandwidth  for  which  the  time  delay  is 
a  constant  value  is  very  small  and  limited.  On  the  other  hand,  in 
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FIG.  5.25  EFFECT  OF  SLAB  THICKNESS  ON  THE  DISPERSION  CURVES. 
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the  region  0  <  t/b  <  0.1,  the  time  delay  approaches  high  values  and 
remains  constant  over  a  large  frequency  range.  This  shows  the  great 
potential  of  thin  films  to  provide  constant  and  high  values  of  group 
time  delay  in  delay  line  applications  with  a  fairly  large  bandwidth. 

Figure  5.27  plots  the  initial  constant  group  time  delay  vs. 
the  normalized  slab  thickness.  In  this  figure,  the  minimum  occurs 
at  t/b  =0.5  and  the  curve  is  unsymmetrical  about  this  minimum  value. 

From  these  observations,  it  can  be  concluded  that  to  obtain 
high  values  of  group  time  delay  over  a  large  bandwidth,  very  thin 
slabs  are  required.  To  increase  the  time  delay  even  more,  it  is 
best  to  choose  a  narrow  width  slab  and  place  it  in  the  center  of  the 
waveguide.  It  should  be  mentioned  that  for  more  accuracy.  Figs. 

5.25,  5>26  and  5.27  are  based  on  the  composition  of  results  obtained 
from  mode  analysis  (Eq.  3.22)  for  large  thicknesses  and  integral 
equation  formulations  (Eqs.  U.t+7)  for  small  thicknesses. 
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FIG.  5.27  INITIAL  CONSTANT  iROUP  TIME  DELAY  VS.  NORMALIZED 


CHAPTER  VI .  CONCLUSIONS 


6.1  Summary  and  Conclusions 

The  objective  of  this  study  vas  to  gain  a  better  understanding 
of  magnetostatic-wave  propagation  in  YIG-loaded  rectangular  wave¬ 
guides.  Both  theoretical  and  numerical  investigations  were  carried 
out  to  achieve  this  objective. 

In  Chapter  II,  the  general  problem  of  wave  propagation  in 
an  unbounded  ferrite  medium  was  formulated  and  magnetostatic  waves 
which  propagate  only  in  a  limited  range  of  wavelengths  were  identi¬ 
fied.  Next,  magnetostatic-wave  propagation  in  bounded  media  was 
studied  and  several  important  structures  were  analyzed.  The  boundary 
conditions  that  these  waves  should  satisfy  at  the  metal  boundaries 
or  on  the  slab  surfaces  were  also  introduced. 

In  Chapter  III,  magnetostatic-wave  propagation  in  a  YIG  slab 
placed  symmetrically  inside  a  waveguide  was  studied.  The  slab  had 
the  same  width  as  the  guide  and  therefore  a  mode  analysis  was  very 
useful  in  formulating  the  dispersion  relations  for  parallel  and 
normal  magnetization.  It  was  also  shown  that  the  obtained  dispersion 
relations  reduced  to  the  published  results  for  the  degenerate  cases 
when  one  or  more  of  the  metallic  boundaries  of  the  guide  were 
driven  to  infinity. 

In  Chapter  IV,  a  modified  configuration  was  studied.  In  this 
configuration  there  was  an  equal  air  gap  on  both  sides  between  the 
slab  and  the  sidewall  of  the  guide.  The  integral  equation  method 
was  effectively  employed  to  obtain  the  dispersion  relations  for  two 
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cases  of  normal  and  parallel  magnetization.  The  obtained  dispersion 
relations  involved  infinite  determinants  which  were  properly  trun¬ 
cated  for  numerical  study. 

In  Chapter  V,  an  extensive  numerical  analysis  was  carried  out. 
Since  most  of  the  obtained  dispersion  relations  were  in  terms  of 
infinite  determinants,  the  proper  truncation  of  these  determinants  was 
studied  and  the  results  of  this  truncation  were  shown.  With  the  help 
of  the  Newton-Raphson  method  the  roots  of  the  truncated  determinant 
were  found  and  plotted  to  obtain  dispersion  and  group  time-delay 
characteristics  for  both  normal  and  parallel  magnetization  cases. 

In  conclusion,  the  analysis  and  numerical  simulation  carried 
out  in  this  work  were  revealing  of  the  magnetostatic-wave  behavior 
in  a  ferrite-loaded  waveguide.  The  width,  thickness,  and  position 
of  the  slab  were  shown  to  be  a  determining  factor  in  the  device 
performance.  The  following  conclusions  can  be  stated  briefly: 

1.  For  parallel  magnetization,  to  obtain  high  time-delay  values, 
the  YIG  slab  must  be  thin,  narrow  and  placed  at  the  bottom  of  the 
guide.  On  the  other  hand,  to  maximize  the  device  bandwidth,  a 
thick,  narrow  YIG  slab  positioned  at  the  top  of  the  guide  is 
preferred. 

2.  For  normal  magnetization,  to  achieve  high  values  of  time 
delay,  again  a  thin,  narrow  slab  must  be  used.  This  time  it  should 
be  placed  in  the  center  of  the  guide.  On  the  other  hand,  high 
bandvidths  can  be  obtained  by  using  wide,  thin  slabs  placed  at  the 
top  or  bottom  of  the  guide. 

It  is  noted  that  there  exists  a  trade  off  between  the:  time 
delay  and  the  device  bandwidth,  an  i  mnximi of.  •'  one  property 
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leads  to  a  poor  value  in  the  other.  Thus  some  design  compromises 
should  be  made. 

In  general,  based  on  the  presented  analysis  in  the  previous 
chapters  and  certain  design  specifications,  an  optimum  device  geometry 
can  be  designed  to  provide  the  best  device  performance  over  a  desired 
frequency  range. 

6.2  Suggestions  for  Further  Study 

The  formulation  developed  in  this  investigation  for  two 
principal  directions  of  magnetization  provides  good  simulation 
results  and  insight  into  the  device  performance.  However,  there  are 
several  additional  topics  that  would  help  to  realize  the  full  poten¬ 
tial  of  magnetostatic-wave  devices  and  thus  need  further  study. 

They  are: 

1.  A  study  of  magnetostatic-wave  propagation  in  a  waveguide  in 
the  presence  of  (a)  a  magnetic  bias  field  in  the  transverse  plane  at 
an  angle  0q  to  the  slab  plane  and  (b)  an  axial  magnetization. 

2.  An  analysis  and  proper  modeling  of  the  propagation  losses. 

3.  An  analysis  of  magnetostatic-wave  propagation  in  a  thick 
YIG  slab  in  a  waveguide  (x  4  0). 

k.  The  design  and  fabrication  of  the  devices  and  experimental 
studies  of  the  magnetostatic-wave  propagation  properties  and 
comparison  with  the  theoretical  results. 

5.  A  study  of  the  higher-order  modes  for  signal  processing 


applications . 


APPENDIX  A.  SAMPLE  COMPUTER  PROGRAMS 

The  following  sample  programs  in  FORTRAN  are  written  to 
solve  for  roots  (K)  of  the  function  G ( K , JMODE ,U ,K1 )  by  using  the 
Nevrton-Raphson  method.  The  function  G(K ,JM0DE,U,K1)  represents 
the  determinant  of  the  truncated  coefficient  matrix.  These  programs 
are  based  on  the  determinant  calculation  procedure  and  the  Newton- 
Raphson  iterative  method  as  described  in  Chapter  V.  The  most 
important  symbols  used  in  this  program  are  defined  as  follows: 

A  Waveguide  width  in  cm. 

B  Waveguide  height  in  cm. 

D  One  half  of  the  slab  thickness  in  cm. 

FRQ  Operating  frequency  in  GHz. 

HI  Internal  dc  magnetic  field  in  Oe. 

JMODE  The  specific  mode  under  consideration. 

K  Wave  number  in  cm-1. 

K1  Off-diagonal  term  of  the  permeability  tensor. 

LN  Size  of  the  truncated  matrix. 

MS  Magnetization  saturation  in  Oe. 

XKX  The  initial  value  of  K  tc  start  the  iteration. 

XO  Gap  length,  in  cm. 

U  Diagonal  term  of  the  permeability  ter.sor. 

DO  Position  of  the  slab  in  the  waveguide  in  cm. 

The  first  samj le  program  finds  the  roots  of  the  dispersion 
relation  as  given  by  the  coefficient  matrix  of  Eq.  where  H , 

is  in  the  x- -direct  ion.  This  program  consists  of  the  foJ lowing  parts: 
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1.  Main  program — finds  the  roots  of  function  G  by  using  the 


Newton-Raphson  method. 

2.  Function  G(K,JM0DE,U,K1 ) — constructs  and  loads  the  truncated 
matrix.  This  function  subprogram  takes  as  its  input  the  wave  number 
(K),  the  specific  mode  under  study  ( JMODE ) ,  and  the  operating  frequency 
(U,K1)  and  returns  the  determinant  of  the  truncated  matrix  in  G  by 
calling  function  XDETtA,N). 

3.  Function  XDET(A,N) — computes  the  determinant  of  the  N  by  N 
truncated  matrix  A  and  returns  the  value  of  the  determinant  in  XDET. 


C . PROGRAM  tt  1 . 

C . PROGRAM  TO  CALCULATE  ROOTS  WHEN  XO  IS  NOT  ZERO  AND 

C . Hdc  IS  PARALLEL  TO  THE  SLAB  PLANE,  i.e.,  Hdc=H0.X 

C . THIS  PROGRAM  CALCULATES  ROOTS  OF  THE  ODD  MODES  ONLY. 

C 

c  **« 

C . THIS  IS  THE  MAIN  PROGRAM  IN  FORTRAN 

C . WRITTEN  TO  SOLVE  FOR  ROOTS  (K>  OF  THE  NONLINEAR 

C . EQUATION  G ( FRQ , K )  =  0  USING  THE  NEWTON-RAPHSON  METHOD, 

C . WHERE  FREQUENCY (FRQ)  IS  A  KNOWN  VALUE. 


IMPLICIT  REAL*4  (K , M) 
DIMENSION  DELTA ( 1 5 ) , K ( 1 5 ) 
COMMON  A , B , D , Z0 , XO , LN 
EPS= , 5 
DK=.  1 
NITER=5 

DELTA ( 1 1=100000. 


KX=- 1 . 

DO  300  NR= 1 ,5 

C . INPUT  STATEMENT  TO  SET  THE  VALUES  OF  FREQUENCY ( FRQ) , 

C . INITIAL  WAVE  NUMBER ( XKX ), GAP  LENGTH ( X0 ), MATRI X  SIZE 

C . ( LN=5 ) , SLAB  POSI TI ON ( Z0 ) , WAVEGUI DE  WI DTH ( A ) , WAVEGUI DE 

C . HEIGHT ( B ) , ONE  HALF  SLAB  THICKNESS ( D ) , SATURATION 

C . MAGNETI ZATION(MS) , INTERNAL  DC  MAGNETIC  FI  ELD (HI )  AND 

C . MODE  UNDER  STUDY  (JMODE )  . 

NAMELI ST/RAD/FRQ , XKX ,  X0 , LN , Z0 , A , B , D , MS , HI , JMODE 
READ  (5, RAD) 

WRITE (6, RAD) 

201  FORMAT ( F 1 0 . 3 , 1  2 ) 

K ( 1 )=XKX 

90  FORMAT ( //3X , 3F20 . 3 ) 

93  FORMAT ( I  2 . 3F 1 0 . 3 ) 

91  FORMAT (2F1 0.3) 


FRQ0'2 .8*HI/1000. 

FRQM=2 . 8*MS/ 1000. 

DEL'FRQO*''  2-FRQ* * 2 
K 1 =FRQ*FRQM/DEL 
U= 1 . +K 1 *FRQ0/FRQ 
WRI TE( 6,90)  U,D,Z0 

C . ROOT  FINDING  BY  ITERATIVE  STEPS 


—  1 3‘>— 


96 


92 

98 

97 


89 

99 
200 

100 
94 
300 


C 

C 

C 


200 


DO  200  1=1, NITER 
11=1+1 

XKDK=K( I ) +DK 

XV=G ( XKDK , JMODE ,  U  ,  K  1  ) 

XW=G ( K ( I ) , JMODE ,  U  ,  K  1  ) 

GPRIM= (XV-XWJ/DK 
WRITE (6,90)  GPRI M , XV , XW 
K ( I  1 ) =K( I ) -XW/GPRI M 
DIFF=K(I 1 ) -K ( I ) 

DELTA(I 1 )=ABS(DIFF) 

XX=G(K(I 1 ) , JMODE, U,K1 ) 

IF  (DELTA (II)  .LE.  EPS)  GO  TO  98 
I F ( DELTA  ill)  .LE.  DELTA ( I  i  )  GO  TO  99 
GO  TO  97 

IF  (ABS(XX)  .LE.  EPS)  GO  TO  100 
GO  TO  92 

WRITE (6, 90)  K ( I ) , DELTA ( I ) , DELTA ( I  1 ) 

K ( 1 ) =K ( 1 ) +KX 
WRI TE ( 6 , 89 )  K ( 1 ) 

FORMAT! 3X, ' K ( 1 )  IS  RESET  TO  A  NEW  VALUE : ' , F9 . 2 ) 

GO  TO  96 

WRITE  (6,94)  I , K ( I  1 ) , DELTA (11), XX 
CONTINUE 

WRITE  (6,94)  I,  K ( I  1 ) , DELTA (II), XX 
FORMAT (I2,2F10.3,F20.3) 

CONTINUE 

STOP 

END 

.FUNCTION  G(F,K)  DESCRIBES  THE  DETERMINANT  OF  AN 
,N  BY  N  TRUNCATED  MATRIX. THE  VALUE  OF  THE  GAP 
,  LENGTH ( X0 )  IS  NOT  ZERO  AND  Hdc  =  H0.X. 

FUNCTION  G ( K , JMODE , U , K 1 ) 

IMPLICIT  REAL ( K ) 

DIMENSION  XM(40,40,2,2),G1(40,2,2),BIG(40,40), BETA (40,40), 
11(40), JJ ( 40 ) 

COMMON  A, B ,D, Z0,X0 , LN 

FORMAT (2X, ' X0= ' ,F7. 3, ’  LN* ' ,12) 

xoa=xo/a 

Z1=Z0-D 
Z2=Z0+D 
PI  =  3  .  14  16 
S  1  =B-Z  1 
S2=B-Z2 

DO  10  I  =  1  , LN , 2 
DO  20  J  = 1  , LN, 2 
ALF  A  =  J  »  PI /A 
ALFA2=ALFA*ALFA 
GAMA  =  SQRT ( ALFA 2  +K  *K ) 

V- 1 . 0 / ( GAMA  *  S I NH ( GAMA  «  5 )  ) 

P2  =  - ( U- 1 . ) * i ALF A 2 ) i ( 2  *  D  +  GAMA*  GAMA ) 

P  1  =  -K 1  *  K - P  2 

P  3  = ( U - ' .  1  *  K  *  K/GAMA 

P4=-P2+K1*K 

C  Z  1  =  COSH  1  GAMA  *  Z 1 ) 

CZ2  =  COSHI GAMA*  Z 2 ) 

CS 1 =COSH ( GAMA*S 1 ) 

CS2--COSH(GAMA*S2) 

G 1 ( J , ’ , I) =V*CZ 1 *( P4+CS 1 *P2*CS2+P3*SINH(GAMA*S 1 ) ) 

G  1  (  J  ,  •  ,-2  )  =V*CZ  1  *  (  P  1  *CS2‘F2*CS1-P3*SINH(  GAMA  *  S2  )  ) 
G1(J,2,i)=V*CS2«(P4*CZ l+P2*CZ2-P3*SINH( GAMA  *  Z 1 ) ) 
G1(J,2,2i=V*CS2»tP'*CZ2*?2*CZ ’  +  P3 *  SI NH ( GAMA  *  Z 2 )  ) 

DO  3  0  L2  =  1  ,  2 
DO  4  0  L  =  1  ,  2 

IF  (I  .EC.  J)  GO  TO  -10 
IJ1=I-J 
I  J  2  =  I  *J 
A1=IJ’*P!»X0A 
A  <:  =  I  J  2  *  P I  *  X  (  A 

BETA  I  I  , J  i =- . b3tb2» (  I S I N ( A  1) /: J 1 )  +  (SIN( A2)/I J2 ) ) 

XM  !  I  ,  J  ,  1.2  ,  L  i  =  BETA  (I  ,  J  )  *G  1  (  J  ,  L2  ,  L  ) 


\J  - 
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50  A3=2*I *PI *X0A 

BETA (I  , X )  =  1 .-2*X0A-SIN( A3 ) / ( I »PI ) 
XM (1,1 , L2 , L ) =  BETA ( 1 , I ) *G 1 ( I , L2 , L ) 
40  CONTINUE 

30  CONTINUE 

20  CONTINUE 

10  CONTINUE 

C . LOADING  THE  BIG  MATRIX 

DO  60  1=1,2 


DO  60  J=  1  ,  2 
DO  55  L= 1 , LN , 2 
II ( L )  =  I +L-  1 
JJ ( L ) = J  +L- 1 

55  CONTINUE 

DO  56  L2= 1 , LN  ,  2 
DO  57  L3= 1 , LN  ,  2 

BIG ( 1 1 ( L2 ) , JJ(L3 ) ) =XM{ L2 , L3 , 1 ,J) 

57  CONTINUE 

56  CONTINUE 

60  CONTINUE 

LN 1 -LN+ 1 
DO  70  1=1 ,LN1 

70  BIG(I , I )  =  BIG(I ,1 )  +  1 .0 

G=XDET( BIG , LN 1 ) 

RETURN 

END 

C . FUNCTION  SUBPROGRAM  TO  CALCULATE  THE  DETERMINANT 

C . OF  AN  N  BY  N  MATRIX  A. 

FUNCTION  XDET ( A , N ) 

DIMENSION  A(40,40) ,B(40,40) 

EPS=. 00001 

C . CONSTRUCT  IDENTITY  MATRIX 

DO  6  I  =  1  ,  N 
DO  5  J=  1  ,N 
IF  (I-J)  4,3,4 

3  B(I,J)=1. 

GO  TO  5 

4  B(I,J)=0.0 

5  CONTINUE 

6  CONTINUE 

C . LOCATE  MAXIMUM  MAGNITUDE  A(I,K)  ON  OR  BELOW  MAIN 

C . DIAGONAL 

DEL= 1 .0 
DO  45  K=  1  , N 
IF  (K-N)  12,30,30 
12  I MAX=K 

AMAX= ABS ( A ( K , K ) ) 

KP 1 =K+ 1 
DO  20  I =KP 1 , N 

IF  ( AMAX-ABS ( A ( I , K ) ) )  15,20,20 
15  IMAX=I 

AMAX= ABS ( A ( I , K ) ) 

20  CONTINUE 

C . INTERCHANGE  ROWS  I  MAX  AND  K  IF  I  MAX  NOT  EQUAL  TO  K 

IF  (IMAX-K)  25,30,25 
25  DO  29  J=  1  ,N 


ATMP  =  A ( I  MAX , J ) 

A ( I MAX , J ) = A ( K , J ) 
A ( K , J ) = ATMP 
BTMP  =  B ( I  MAX , J ) 

B ( I  MAX , J ) =  B ( K , J ) 


29  B ( K , J )  =  BTMP 
DEL=-DEL 

30  CONTINUE 

IF  (ABS(A(K,K) ) -  EPS )  93,93,35 

C . TEST  FOR  SINGULAR  MATRIX 

35  DEL= A ( K , K ) *DEL 

C . DIVIDE  PIVOT  ROW  BY  ITS  MAIN  DIAGONAL  ELEMENT 

DI V= A ( K , K ) 

DO  38  J=  1  ,N 
A  (  K  ,  J  )  =  A  (  K  ,  J  )  /D I V 
38  B(K,J)=B(K,J)/DIV 

c . REPLACE  EACH  ROW  BY  LINEAR  COMBINATION  WITH  PIVOT  ROW 
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39 

DO  43  1  =  1  ,N 

AMULT= A ( I , K ) 

IF(I-K)  39,43,39 

DO  42  J= 1 ,N 

42 

A  ( I  ,  J  )  =  A  ( I  ,  J  ) -AMULT*A ( K , J ) 
B ( I , J)=B(I , J)-AMULT*B(K, J) 

43 

CONTINUE 

45 

CONTINUE 

99 

XDET=DEL 

RETURN 

93 

GO  TO  99 

END 

The  second  sample  program  finds  the  roots  of  the  dispersion 
relation  for  zeroth-  and  first-order  modes  as  described  by  a  5  by  5 
matrix  given  by  Eq.  5.10  where  is  in  the  z-direction.  This 
program  has  the  following  parts: 

1.  Main  program — finds  the  roots  of  function  G. 

2.  Function  0(K, JM0DE,U,K1 ) — constructs  and  loads  the 
truncated  5  by  5  matrix.  This  function  subprogram  takes  as  its 
input  the  wave  number  (K),  the  specific  mode  under  study  (JMODE),  and 
the  operating  frequency  (U,Kl)  and  returns  the  determinant  of  the 
truncated  matrix  in  G  by  calling  subroutine  XIUV(A,B,XDET,N) . 

3.  Subroutine  XINV(A,B,XDET,N) — computes  the  inverse  of  the 
u  by  bi  truncated  matrix  A  and  returns,  it  in  matrix  B.  This 
subroutine  also  computes  the  determinant  of  mati ix  A  and  returns 
it  in  XDE’T. 


C . PROGRAM  #2 . 

C . PROGRAM  TO  CALCULATE  ROOTS  WHEN  XO  IS  NOT  ZERO  AND 

C . Hdc  IS  NORMAL  TO  THE  SLAB , i . e . ,  Hdc=H0.Z 

C . THIS  PROGRAM  CALCULATES  ROOTS  OF  THE  FIRST  MODE  ONLY. 

C  *** 

C . THIS  IS  THE  MAIN  PROGRAM  WRITTEN  IN  FORTRAN  SOLVING 

C . FCR  ROOTS  (K)  OF  THE  NONLINEAR  EQUATION  G(FRQ,K)=0. 


C . USING  THE  NEWTON- P.APHSON  METHOD,  WHERE  FRQ  IS  KNOWN. 

IMPLICIT  REAL*4  (K,M) 

DIMENSION  DELTA (  I  5 ) , K (  ! 5 ) 

COMMON  A,B,D,Z0,X0,LN,MS,HI 
EPS- . 5 
DK=  .  1 
NITERS 

DELTA ( 1 )= 100000. 

KX=- 1 . 

no  nn  hr=  i  ,  5 
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. INPUT  STATEMENT  TO  SET  THE  VALUES  OF  FREQUENCY ( FRQ) , 
.INITIAL  WAVE  NUMBER < XKX ), GAP  LENGTH ( XO ), MATRIX  SIZE 
. ( LN  =  5 )  , SLAB  POSITION ( ZO ) , WAVEGUI DE  WI DTH ( A ) .WAVEGUIDE 
.HEIGHT(B) ,ONE  HALF  SLAB  THICKNESS (D) .SATURATION 
. MAGNETIZATION (MS) .INTERNAL  DC  MAGNETIC  FI ELD(HI ) .MODE 
.UNDER  STUDY ( JMODE ) . 

NAMELI ST/RAD/FRQ , XKX .  XO , LN , ZO , A , B , D .MS , HI , JMODE 
READ  (5, RAD) 

WRITE ( 6 , RAD ) 

FORMAT ( F 1 0 . 3 , 1  2 ) 

K ( 1 )=XKX 

FORMAT  (//3X.3F20.3) 

FORMAT ( I  2 , 3F 1 0 . 3 ) 

FORMAT (2F1 0.3) 

FRQ0=2.8*HI/1000. 

FRQM=2 . 8»MS/1 000 . 

DEL=FRQ0»*2-FRQ**2 
K 1 =FRQ*FRQM/DEL 
U=1 . +K 1 *FRQ0/FRQ 
WRI TE ( 6 , 90 )  U.K1.Z0 
.ROOT  FINDING  BY  ITERATIVE  STEPS 
DO  200  1=1 .NITER 
11=1+1 

XKDK=K ( I ) +DK 

XV=G  ( XKDK  ,  JMODE  ,  U  ,  K  1  ) 

XW=G ( K ( I ) , JMODE, U, K1  ) 

GPRIM=(XV-XW)/DK 
WRITE (6,90)  GPRIM.XV.XW 
K ( I  1 ) =K ( I ) -XW/GPRI M 
DI FF=K ( I  1 ) -K { I ) 

DELTA(I 1 ) = ABS ( DI FF ) 

XX=G ( K ( I  1 ) , JMODE, U,K1 ) 

IF  ( DELTA (II)  .LE.  EPS)  GO  TO  98 
IF (DELTA ( I  1 )  .LE.  DELTA ( I ) )  GO  TO  99 
GC  TO  97 

IF  ( ABS ( XX )  .LE.  EPS)  GO  TO  100 
GO  TO  92 

WRITE (6,90)  K ( I ) ,DELTA(I ) .DELTA (I  1 ) 

K  (  1 ) *K( 1 ) +KX 
WRITE( 6,89)  K (  1  ) 

FORMAT ( 3X , ' K ( 1 )  IS  RESET  TO  A  NEW  VALUE : ' , F9 . 2 ) 

GO  TO  96 

WRITE  (6,94)  I ,K(I 1 ) ,DELTA(I 1 ) ,XX 
CONTINUE 

WRITE  (6,94)  I,  K(I 1 ) ,DELTA(I 1) ,XX 
FORMAT (I 2 , 2F 1 0 . 3 , F20 . 3 ) 

CONTINUE 

STOP 

END 

.FUNCTION  G(F,K)  DESCRIBES  THE  DETERMINANT  OF  A 
.5  BY  5  MATRIX  FOR  THE  FIRST-ORDER  MODE,  WHEN 
.Hdc=H0.Z  AND  X0  IS  NOT  EQUAL  TO  ZERO. 

FUNCTION  G (K, JMODE, U.K1  ) 

IMPLICIT  REAL ( K ) 

DIMENSION  XA (40,40) ,XB(40,40) ,VN(2) , UN( 2 ) , BIG ( 40 , 40 ) 

COMMON  A,B,D,Z0,X0,LN, MS , HI 

Z 1=Z0-D 

Z2=Z0+D 

PI =3.  14  16 

BZ0=B-Z0 

BZ 1  =  B- Z 1 

BZ2=B-Z2 

X0A=X0/A 

AX0=A-X0 

A2X0=A-2*X0 

A0=-PI*X0/A2X0 

A1=-2*A0/A 

A2«2*PI « ( 1-X0AJ/A2X0 
PA2=PI *AX0/A 
PA  1 =PI *X0A 
A3=A0+A2*AX0 


n  n 
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A4=A0*A2*X0 
A5=A0*A1 *AXO 
A6=A0+A1*X0 

S34=.5*(SIN(A3)-SIN(A4) )/A2 
S56-.5*(S:N(A5)-S:N(A6) )/A' 

C34= . 5* 1 -COS ( A3 ) +COS ( A4 ) )/A2 
C56= . 5* ( -COS ( A5 ) +COS ( A6 ) ) /A  1 
XA (1,1 ) =  A2X0 
XA( 1 , 2)=0. 

XA ( 1 , 3)=2»A2X0/PI 

XA ( 2 ,  1  )  =  . 3  1 8* A* ( SI N ( PA2 ) -SIN ( PA  1 ) ) 

XA(2,2)=S34*S56 
XA (  2 , 3 ) =C 34  +C 56 

XA (3,1)=. 318*A*(COS(PA2) -COS ( PA  1 ) ) 

XA(3,2)=C34-CS6 
XA(3,3)=S34~S56 
NN=  3 

CALL  XINV(XA,XB,XDET,NN) 

DO  10  1=1,2 
11=1-1 

ALFA=II*PI/A 
ALFA2= ALFA* ALFA 
GAMA  =  SQRT (ALFA2+K*K) 

W=1 . 0/ ( GAMA* S I NH (GAMA*  B ) } 

F=(U-1 . )/(U*D*D) 

SZ0  =  SINH( GAMA*  Z0 ) 

SZ 1=SINH(GAMA*Z 1 ) 

SZ2*SINH ( GAMA*Z2 ) 

SBZ0=SI NH ( GAMA*BZ0 ) 

SBZ 1 =  SINH (GAMA*BZ 1 ) 

SBZ2=SI NH ( GAMA*BZ2 ) 

CZO=COSH ( GAMA*Z0 ) 

CZ  1=COSH(GAMA*Z1 ) 

CBZO=COSH (GAMA*  BZO } 

CBZ2=COSH(GAMA*BZ2) 

B1=(SZ0-SZ1 ) /GAMA 
B2=(SZ2~SZ0) /GAMA 
B3=(-SBZ0+SBZ D/GAMA 
B4-(-SBZ2*SBZ0) /GAMA 

VN ( I ) =  W* ( ( B 1 *B2) *CBZ2+ (B3+B4)*CZ1“2*B1 *CBZ0- 2*B4 *CZ0 ) 
UN (I ) =  VN ( I ) +4*W*B 1 *CBZO  +  4*W*B4*CZO 
10  CONTINUE 

CC1=COS(PA 1 ) 

CC2=COS ( PA2 ) 

DC0=-.25*(K*K*K1*PI  )*(CC1+CC2)/A2X0 
U 1 4= ( U" 1 )/4 

DC  UU14* (XB<  1  ,  1 ) +XB ( 2  ,  1  )  ) 

DC2  =  U14* (XB(  1 ,  1 ) -XB  <  2 ,  1  )  ) 

DC3=U 14* (XB( 1 , 2) +XB(2, 2) ) 

DC4  =  U14*(XB< 1 ,2)-XB(2,2)  ) 

DC5  =  U 1 4  * { XB ( 1 ,  3)+XB(2, 3) ) 

DC6=U 1 4  * ( XB (  1,3) -XB (2,3) ) 

D0=-DC 1 *CC 1 +DC2*CC2+DC0*XB ( 3 , 1) 

D1  =  -DC3*CC1 *DC4  *CC2+DC0  *XB (3,2) 

D2=-DC5*CC1 +DC6 *CC2*DC0 *XB  (  3 , 3 ) 

A00= 1 . -2*X0A 

A0l=.636b*(SIN(PA2)-SIN(PA1) ) 

A  10  =  AO  1/2. 

A1 1=A00+.  1 59*(SIN( 2*PA2)-SIN( 2*PA1  )  ) 

E10=.316*(COS(PA1 ) -COS ( PA2 ) i 
B 1  1  =  .  159* ( COS (2* PA  1 ) -COS ( 2  *  PA2 ) ) 

AVO  =  AO  1  * VN ( 2  ) 

AVI =A 1  1  *  VN ( 2 ) 

AU0  =  AG  1 *UN ( 2 ) 

AU 1 = A  1  1*UN(2) 

BU 1 =  B 1  1 *UN ( 2  ) 

LOADING  THE  MATRIX  IN  THE  FINAL  FORM  READY 
FOR  DETERMINANT  CALCULATION. 

B I G (  1,1)  =  ! . +  F  *  A  0  0  *  VN (  1 ) 

BIG (  1 ,2)=F*AV0 
B I G (  1,3) =D0*AV0 
B I G (  1 , 4 ) =  D 1  * AVO 
BIG( 1,5) =D2  * AVO 
B I G  (  2  ,  1  )  =  F  *  A  1  0  *  VN  (  1  ) 


o  fo  -*  —  onff>ui>t*  w  n  o  o  non 
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BIG ( 2 , 2 ) = 1 .+F*AV1 
BIG( 2 , 3 ) »D0«AV1 
BIG ( 2 , 4 ) =D 1 *AV 1 
BIG( 2 , 5 ) =D2«AV1 
BIG( 3 , 1 )=F*AOO*UN(  1) 

EIG ( 3 , 2 ) =F*AU0 
BIG ( 3 , 3 ) = 1 . +D0»AU0 
BIG ( 3 , 4 ) *D 1 *AU0 
BIG( 3,5) =D2»AU0 
BIG( 4 ,  1 ) »F*A 1 0*UN (  1} 

BIG(4,2)-F*AU1 
BIG ( 4 , 3 ) *D0*AU1 
BIG( 4 , 4 ) = 1 . +D1*AU1 
BIG(4,5)=D2»AU1 
BIG ( 5 , 1 )=F»B10*UN(  1 ) 

BIG( 5 , 2 ) =F*BU1 
BIG ( 5 , 3 ) =D0*BU 1 
BIG ( 5 , 4 ) =D1 *BU1 
BIG{ 5 , 5) = 1 .+D2*BU1 
LN=5 

CALL  XI NV ( BIG , XB , XDET , LN ) 

G=XDET 

RETURN 

END 

. SUBROUTINE  FOR  DETERMINANT  CALCULATION  AND 

. MATRIX  INVERSION  BY  ELIMINATION  WITH 

. PARTIAL  PIVOTING. 

SUBROUTINE  XINV< A , B , XDET, N ) 

. ABORIGINAL  MATRIX,  B=IVERSE  MATRIX,  XDET= 

. DETERMINANT,  N=MATRIX  SIZE. 

DIMENSION  A(40,40) , B ( 40 , 40 ) 

EPS=. 00001 

. CONSTRUCT  IDENTITY  MATRIX  B ( I , J )  =  1 

DO  6  I «  1  ,  N 
DO  5  J-1.N 
IF  (I-J)  4,3,4 
B(I ,J)=1 . 

GO  TO  5 
B ( I , J ) =0 . 0 
CONTINUE 
CONTINUE 

. LOCATE  MAXIMUM  MAGNITUDE  A(I,K)  ON  OR  BELOW  MAIN 

. DIAGONAL 

DEL= 1.0 
DO  45  K=  1  , N 
IF  (K-N)  12,30,30 
2  IMAX=K 

AMAX= ABS ( A ( K , K )  ) 

KP1=K+1 
DO  20  I =KP 1 , N 

IF  ( AMAX- ABS ( A ( I , R ) ) )  15,20,20 
5  I  MAX* I 

AMAX= ABS ( A ( I , K ) ) 

0  CONTINUE 

. INTERCHANGE  ROWS  IMAX  AND  K  IF  I  MAX  NOT  EQUAL  TO  K 

IF  (IMAX-K)  25,30,25 

25  DO  29  1 ,N 

ATMP=A ( IMAX , J ) 

A ( I  MAX , J ) = A ( K , J ) 

A(K,J)=ATMP 
BTMP=  B ( I  MAX , J ) 

B ( I MAX , J ) =  B ( K , J ) 

29  BIK,J)*BTMP 
DEL=-DEL 

30  CONTINUE 
IF  (ABS ( A ( K , K ) ) -EPS )  93,93,35 
TEST  FOR  SINGULAR  MATRIX 
DEL=A( K , K ) *DEL 


C. 

35 


V 
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c 


3S 

C. 


39 

42 

43 
45 

99 

93 


DIVIDE  PIVOT  ROW  BY  ITS  MAIN  DIAGONAL  ELEMENT 
DI V= A ( K , K ) 

DO  38  J=1,N 
A(K,J)=A(K,J)/DIV 

.replace^ach'row  by  LINEAR  COMBINATION  WITH  PIVOT  ROW 
DO  43  1=1, N 
AMULT=A( I ,K) 

IF(I-K)  39,43,39 
DO  42  J=1,N 

A(I,J)-A(I,J) -AMULT* A ( K , J ) 

TUT  ,J)=B(I . J)-AMULT*B(K,J> 


CONTINUE 

CONTINUE 

XDET=DEL 

RETURN 

GO  TO  99 

END 
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